Astronomy & Astrophysics manuscript no. paper 


©ESO 2012 


March 2, 2012 





The Galactic WC stars 
Stellar parameters from spectral analyses indicate a new evolutionary sequence 

A. Sander, W.-R. Hamann, and H. Todt 

Institut fur Physik und Astronomie, Universitat Potsdam, Karl-Liebknecht-Str. 24/25, D-14476 Potsdam, Germany 
e-mail: ansander@astro . physik . uni -po t sdam . de , wrh@astro . physik . uni -potsdam . de 

Received 5 August 201 1; accepted 30 January, 2012 

ABSTRACT 

Context. The life cycles of massive stars from the main sequence to their explosion as supernovae or gamma ray bursts are not 
yet fully clear, and the empirical results from spectral analyses are partly in conflict with current evolutionary models. The spectral 
analysis of Wolf-Rayet stars requires the detailed modeling of expanding stellar atmospheres in non-LTE. The Galactic WN stars have 
been comprehensively analyzed with such models of the latest stage of sophistication, while a similarly comprehensive study of the 
Galactic WC sample remains undone. 

Aims. We aim to establish the stellar parameters and mass-loss rates of the Galactic WC stars. These data provide the empirical basis 
of studies of (i) the role of WC stars in the evolution of massive stars, (ii) the wind-driving mechanisms, and (iii) the feedback of WC 
stars as input to models of the chemical and dynamical evolution of galaxies. 

Methods. We analyze the nearly complete sample of un-obscured Galactic WC stars, using optical spectra as well as ultraviolet spectra 
when available. The observations are fitted with theoretical spectra, using the Potsdam Wolf-Rayet (PoWR) model atmosphere code. 
A large grid of line-blanked models has been established for the range of WC subtypes WC4 - WC8, and smaller grids for the WC9 
parameter domain. Both WO stars and WN/WC transit types are also analyzed using special models. 

Results. Stellar and atmospheric parameters are derived for more than 50 Galactic WC and two WO stars, covering almost the whole 
Galactic WC population as far as the stars are single, and un-obscured in the visual. In the Hertzsprung-Russell diagram, the WC stars 
reside between the hydrogen and the helium zero-age main sequences, having luminosities L from 10 4 - 9 to 10 5 6 L B . The mass-loss 
rates scale very tightly with L 8 . The two WO stars in our sample turn out to be outstandingly hot («;200kK) and do not fit into the 
WC scheme. 

Conclusions. By comparing the empirical WC positions in the Hertzsprung-Russell diagram with evolutionary models, and from 
recent supernova statistics, we conclude that WC stars have evolved from initial masses between 20 solar masses and 45 M Q . In 
contrast to previous assumptions, it seems that WC stars in general do not descend from the most massive stars. Only the WO stars 
might stem from progenitors that have been initially more massive than 45 M . 

Key words. Stars: evolution - Stars: mass-loss - Stars: winds, outflows - Stars: Wolf-Rayet - Stars: atmospheres - Stars: massive 



1. Introduction 

Wolf-Rayet (WR) stars represent advanced evolutionary stages 
in the evolution of massive stars. Owing to their high tem- 
peratures and high mass-loss rates, WR stars provide a large 
amount of ionizing photons, mechanical momentum, and mat- 
ter to their environment, making WR stars a keystone in the cos- 
mic evolution. Because of the extreme non-local thermal equilib- 
rium (non-LTE) situations in the expanding atmospheres, their 
spectral analysis only became possible after adequate, power- 
ful and complex model atmosphere codes had been established. 
Since the early 90s, these tools ha ve been available. Early anal- 
yses h a ve been published, e.g ., bv lHillieJ dl987h.lSchmutz et alJ 
Jl989l) . lHamann et all d 19931) . and lHamann et alJ d 1995b for the 
WN stars, i.e. WR stars with strong nitrogen emission. For 
the WC subclass, which are characterized by strong carbon 
emissi o n-lines, wind models have been applied by, e.g., iHillierl 
( 19891. lHamann et alJ dl992h . and lKoesterke et alJdl992l) . 

The models have since been significantly improved. With 
the inclusion of iron line blanketing and microclumping 
dHillier & Milled Il999t iGrafener et all 120021) . a new genera- 
tion of models have allowed better spectral fits, and led to 
a revision of the empirical stellar parameters. A comprehen- 



sive re-analys is of th e Galactic WN stars was performed by 
Ha mann et alJ d2006l) . Basically all W N stars in the Large 
Magel lanic Cloud have been analyzed by Ruhling (2008, and in 
prep.l. lPasemann et alJ d201 ll) studied the WR stars in the Small 
Magellanic Cloud. In contrast to the WN spectra, the WC spec- 
tra turned out to be much harder to reproduce. Our model spectra 
can fit most, but not all of the observed features, leaving some 
ambiguity about the choice of the best-fit model. 

The present study is based on models that are not hydro- 
dynamically consistent. The mass-loss rate and wind veloc- 
ity field are instead treated as free parameters. The only hy- 
drodynamical WC model that has been constructed so far is 
that developed to study th e prototypical WC5 star WR 111 
( Graf ener & Hamannl [2005). This model arrived at a much 
higher stellar temperature (140kK) and smaller mass-loss rate 
than the semi -empirical approach o f the present paper. Using this 
same model, Grafener et d] ((2002) had obtained T, = 85 kK for 
WR 1 1 1 , while only 62 kK had been derived with non-blanketed 
models dKoesterke & Ha mann 1995]). 

In the afore mentioned paper of iKoesterke & Hamannl 

d 19951) . 25 Galactic WC stars of subtypes WC4 to WC8 were 
analyzed with models containing only carbon and helium. After 
the inclusion of oxygen and the line-blanketing by iron-group 



1 



A. Sander et al.: The Galactic WC stars 



elements, Barniske et al.l (120061) re-analyzed a few of these WC 
stars and found that a carbon mass fraction of 40% appears to be 
typical. 

The present study now covers for the first time a large part 
of the Galactic sample of WC stars of all subtypes. The analy- 
ses comprise 40 WC stars of subtypes WC4 to WC9. In addition, 
we consider six stars of the WN/WC transition type and two WO 
stars. The spectral analyses presented in this work are appropri- 
ate for single-star spectra. A few binaries are included in order to 
demonstrate the bias of the obtained parameters if the composite 
nature of their spectra is neglected. 

The empirical stellar parameters and mass-loss rates of the 
Galactic WC stars established by our analyses provide the em- 
pirical basis for studies of (i) the role of WC stars in the evolution 
of massive stars, (ii) the wind-driving mechanisms, and (iii) the 
feedback of WC stars as input to the modeling of the chemical 
and dynamical evolution of galaxies. 

In the next section, we introduce our sample of stars and the 
observational data used in this work. The model calculations are 
characterized in Sect.[3] In Sect.|H we describe the analyses, and 
in Sect.|5]we present the results and comment on individual stars. 
Sect.|6]then discusses the results with regard to stellar evolution, 
before arriving at our conclusions in the final section (Sect.|7]i. 

2. Observational data 

The current WR catalog from Ivan der Huchil(l200"li r2006) lists 
298 WR stars in our Galaxy, of which 113 are of WC type. 
However, most of the WR sta rs that have been add ed to the cata- 
log since its previous release (Ivan der H ucht 2001) were discov- 
ered in infrared observations, and suffer from high interstellar 
extinction in the visual range. Since our analyses are based on 
optical spectra, our sample is basically restricted to the list of 
Ivan derHuchid200lh . 

A significant fraction of the WR stars are known or sus- 
pected to have a binary companion, mostly of O typ e. Some 
of the stars analyzed by Koeste rke & Hamannl (1 19951) . for in- 
stance, have since been found to be binaries. In WR+O bina- 
ries, the spectrum is composed of both components. Often the 
O-star's contribution even dominates the overall spectrum. The 
decomposition of these spectra is beyond the scope of this pa- 
per. Nevertheless, we include a few of these composite spectra 
in our study and analyze them as if they were single WR stars, 
to illustrate the effects on the results. 

A visible effect of an O-star companion is the so-called "di- 
luted emission lines" (d.e.l.). When the light of a companion 
contributes significantly to the spectral continuum, the emission 
lines appear weaker in the normalized spectrum, hence are called 
"diluted". However, we do not a priori judge a WC star as a 
binary on the basis of apparently diluted emission lines alone. 
The proof of binarity can only be derived from radial-velocity 
variations or an eclipse light curve. The d.e.l. binary criterion is 
based on the assumption that stars of the same su btype have sim- 
ilarly strong emission lines. For WN stars. lHamann et al. I d2006h 
demonstrated that this is not the case. However, we see below 
that the WC sequence has rather uniform spectra, making the 
d.e.l. criterion a valid binary indicator. 

For a few putatively single, visually un-obscured WC stars, 
we have no spectra at our disposal, which leaves a sample of 
56 stars for our study. In addition to the proper WC stars, our 
sample comprises two WO stars, five WN/WC transition types, 
and one WC/WN star (see TableQ}- 

For the analyses, w e collected optical spectra from both 
iTorres & Massevl (U987) and our own observations performed 



at the Calar Alto observatory or ESO. The optical spectra typ- 
ically cover the wavelength range 3300-7300 A. In some cases 
however, if a star e.g. is not bright enough in a certain wave- 
length regime or there were technical difficulties during the ob- 
servation, only a part of this range is covered. The typical spec- 
tral resolution is around 10 A (ITorres & Massevlll987l) . For the 
ultraviolet range, IUE spectra were retrieved from the INES 
archive. They cover the range from 11 50 A to 3200 A with a 
typical resolution of around 6 A. Inf rared photometry (J, H , K) 
was taken from the 2MASS catalog (ISkrutskie et al.ll2006l) . and 
in two cases (WR52, WR 142) was augmented by Spitzer IRS 
observations. For so me stars, we used additional opti cal narrow- 
band photometry b y iLundstrom & Stenholml ([1984), MSX in- 
frare d photometry dEgan et al.l 120031) . and/or Spitzer spectra 
from lArdila eTal] ( 1201 Ol) . 

3. Model grid 

To reproduce the observed spectra, we used the Potsdam 
Wolf-Rayet (PoWR) code for expa nding stellar a t mosph eres. 
The model basics are desc ribed in iGrafener et al.l d2002l) and 
lHamann & Grafenerl d2004l) . The non-LTE radiative transfer is 
calculated for a spherically symmetric and stationary outflow 
with a pre-described velocity field. In the supersonic part, we 
adopt the usual /3-law with the terminal velocity Voo- The expo- 
nent /3 is set to 1 .0 throughout this study, except for the two WO 
stars for which /3 = 0.8 gave slightly better fits. In the subsonic 
part, the velocity field is defined by the requirement that the den- 
sit y stratification approaches the hydrostatic limit, as described 
bylHamann & Sch mutzl d!987l) . Since the mass-loss rate M and 
the velocity field are treated as free parameters, the models are 
not hydrodynamically consistent. 

Wind inhomogeneities are taken into account in the "micro- 
clumping" approximation, assuming that individual "clumps" 
are optically thin at all frequencies, and that the interclump 
medium is void. The clumping factor D denotes the density 
within the clumps, relative to the homogeneous model with the 
same M, and is assumed to be constant throughout the atmo- 
sphere in this study. 

Complex model atoms are taken into account, as summarized 
in Table|2] Ir on-group elements ar e treated in the superlevel ap- 
proximation (Grafener et al. 2002). 

The inner boundary of our model atmospheres is located at a 
radial optical depth of tr oss = 20, which we define to correspond 
to the stellar radius R*. We define the "stellar temperature" T, 
as the effective temperature related to that radius, via Stefan- 
Boltzmann's law 

L = 4nRlo-sBT* (1) 

and th e stellar lumino s ity L. 

As ISchmutz et al.l ([1989) found out, stellar winds produce 
very similar emission-line spectra when they agree in a certain 
combination of the mass-loss rate and other parameters. To ex- 
ploit this scaling invariance, they introduced the so-called "trans- 
formed radius" 



1 *[2500km/s/ 10- 4 M o /yrj 

(th e clumping factor D has be en incorporated in this definition 
by lHamann & Koes terke 199a). 

The name "tran sformed radius", historically coined by 
ISchmutz et al.l dl989l) . is actually misleading since R t does not 
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Table 1. Galactic WC stars, not visually obscured (Ivan der Huchtll2001h 



WR 


Subtype 1 


Figure 2 


Remarks 3 


WR 


Subtype 1 


Figure 2 


Remarks 3 


4 


WC5 


B.l 












5 


WC6 


B.2 




75a 


WC9 


- 


no spectra 


7a 


WN4h/WC 


- 


no spectra 


76 


WC9d 


- 


no spectra 


8 


WN7/WCE 


B.3 


binary WN+WC? 


77 


WC8+OB 


- 


binary 


9 


WC5+07 


- 


binary 


79 


WC7+05-8 


- 


binary 


11 


WC8+07.5 


- 


binary 


80 


WC9d 


B.27 




13 


WC7 


B.4 




81 


WC9 


B.28 




14 


WC7+? 


B.5 


binary? no d.e.l. 


86 


WC7+B0I-III 


B.29 


binary: pseudo fit 5 


15 


WC6 


B.6 




88 


WC9+? 


B.30 


WN+WC or WN/WC9 6 


17 


WC5 


B.7 




90 


WC7 


B.31 




19 


WC4pd+09.6 


- 


binary 


92 


WC9 


B.32 




20a 


WN7:h/WC 


- 


now WN6ha+WN6ha 4 


93 


WC7+07-9 


- 


binary 


23 


WC6 


B.8 




95 


WC9d 


B.33 




26 


WN7/WCE 


B.9 




96 


WC9d 


- 


no spectra 


27 


WC6 


B.10 




98 


WN8/WC7 


B.34 


now WN7/WC+08-9 7 : pseudo fit 


30 


WC6+06-8 


- 


binary 


98a 


WC8-9vd+? 


- 


no spectra 


30a 


W04+05-5.5 


- 


binary 


101 


WC8 


- 


no spectra 


31c 


WC6+OB 


- 


binary 


102 


W02 


B.35 




32 


WC5+OB 


- 


binary 


103 


WC9d 


B.36 




33 


WC5 


B.ll 




104 


WC9d+B0.5V 


B.37 


binary: pseudo fit 5 


38 


WC4 


B.12 


binary? 


106 


WC9d 


B.38 




38b 


WC7+OB 


- 


binary 


107a 


WC6 


- 


no spectra 


39 


WC7+OB 


B.13 


binary: pseudo fit 5 


111 


WC5 


B.39 




41 


WC5+OB 


- 


binary 


112 


WC9d+OB 


- 


binary 


42 


WC7+07V 


- 


binary 


113 


WC8d+08-9IV 


B.40 


binary: pseudo fit 5 


45 


WC6 


B.14 




114 


WC5 


B.41 




47c 


WC5 


- 


no spectra 


117 


WC9d 


B.42 




48 


WC6+09.5 


- 


binary 


118 


WC9d 


- 


no spectra 


48a 


WC9d 


- 


no spectra 


119 


WC9d 


B.43 




48b 


WC8ed 


- 


no spectra 


121 


WC9d 


B.44 




50 


WC7+OB 


B.15 


binary: pseudo fit 5 


125 


WC7ed+09III 


B.45 


binary: pseudo fit 5 


52 


WC4 


B.16 




126 


WC5/WN 


B.46 




53 


WC8d 


B.17 




132 


WC6 


B.47 




56 


WC7 


B.18 




135 


WC8 


B.48 




57 


WC8 


B.19 




137 


WC7pd+09 


B.49 


binary: pseudo fit 5 


58 


WN4/WCE 


B.20 




140 


WC7pd+04-5 


- 


binary 


59 


WC9d+OB? 


B.21 


binary? 


142 


W02 


B.50 




60 


WC8 


B.22 




143 


WC4+Be 


B.51 


binary: pseudo fit 5 


64 


WC7 


B.23 




144 


WC4 


B.52 




65 


WC9d+OB? 


B.24 


binary? 


145 


WN7/WCE 


B.53 




68 


WC7 


B.25 




146 


WC6+08 


B.54 


binary: pseudo fit 5 


69 


WC9d+OB? 


B.26 


binary? 


150 


WC5 


B.55 




70 


WC9vd+B0I 




binary 


153ab 


WN6/WCE+O6I 




binary 


73 


WC9d 




no spectra 


154 


WC6 


B.56 





Notes. {1) Classification from van der Hucht (2001) (2) Spectral fits, shown in the online appendix (3) Further comments to the individual stars 
ar e given in Sect.[Al ( 4 ) iNaze et alj d2008h (5) Known binary fitted as if it were a single star (as discussed in Sect.[2j (6) Will iams et alJ d2005h 
(7) iGamen & Niemelal d2003fi 



have the meaning (although the units) of a radius. A more appro- 
priate quantity to consider would be R~ 3 , which might be called 
a "normalized emission measure". Since this quantity is propor- 
tional to the volume integral of the density squared, divided by 
the stellar surface, Z?J~ 3 scales with the emission from recombina- 
tion lines normalized to the continuum. This explains why dif- 
ferent combinations of R*, »„, and mass-loss rate M result in 
approximately the same normalized WR emission-line strengths 
as long as R t (or R~ 3 ) is kept at the same value. 

Wolf-Rayet model spectra are most sensitive to T* and R t . 
When finding systematically the model that provides the clos- 
est fit to a given observation, it is therefore most convenient to 
rely on model grids in which these two parameters are varied. 
The other parameters, namely the chemical composition and the 
terminal wind velocity, are kept constant within one model grid. 



We established one large model grid for most of the WC pa- 
rameter space, and two smaller grids for the WC9 domain. The 
grid resolution is 0.05 in log T* and 0. 1 in log R t . For the chem- 
ical composition, we assumed mass fractions of 55% helium, 
40% carbon, 5% oxygen, and 0.16% iron-group elements (cf. 
Table|2]l. To verify that these abundances are adequate to the ob- 
served stars, we also calculated smaller grids with 60% and 20% 
carbon fractions, respectively. 

For the main WC grid, we set the terminal wind velocity to 
2000 km s~', which is in the middle of the range actually ob- 
served. The two WC9 grids are calculated for lOOOkms 1 and 
1600kms~'. Owing to their their lower ionization, the WC9 
models include an extended C n model atom. 

Throughout this study, the microturbulence velocity is set to 
lOOkm/s, based on the experience that this yields profile shapes 
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similar to the observed ones. For the clumping density factor, 
we adopt D — 10, which generally reproduces the observed 
strength of the electron-scattering wings. However, the true de- 
gree iof_tiie_W2ndinhomogeneity remains the subject of debate 
(cf. lHamann et al.l 12008). According to Eq. ©, the empirically 
obtained mass-loss rates scale with the adopted clumping factor 
as M oc ZT 1/2 . 

The luminosity was kept fixed over all our grids at 
log L/L Q = 5.3. According to Eqs. ([TJ and (fj), the scaling of the 
model to a different luminosity implies a scaling of the mass-loss 
rate as M oc l}l A . In a strict sense, the scaling invariance between 
models of the same R t is only approximate, but accurate enough 
to cover the range of luminosities in our sample. 

Figures Q] and [2] display contour plots for the strength of the 
brightest WC emission lines. The small dots represent the cal- 
culated models of the standard WC grid. Each contour refers to 
a certain equivalent width, labeled with —Wa in A. The ratio of 
the lines Cm 5696 A to Civ 5808 A is u sed as the primary cri - 
terion for the WC subtype classification dCrowther et al.lfl998h . 
For ea rly subtypes, O v 5590 A is also considered (Ivan der Huchtl 
12001 . 

For the first time, we apply the PoWR code when analyz- 
ing spectra of the late subtype WC9. Most of these stars show 
an excess of infrared emission that is attributed to circumstellar 
dust, as denoted by the letter "d" appended to their classification, 
WC9d. This circumstellar dust usually leads to strong reddening 
of the stellar radiation, inhibiting their observation in the UV. 
The optical spectrum, however, does not suffer much from the 
dust emission and absorption, and can be analyzed with the help 
of our models together with the non-dusty WC9 stars. 

The WC9 stars have a lower stellar temperature than the ear- 
lier WC subtypes, which is evident from the lower ions that are 
visible in their spectra. We therefore calculated special, small 
grids of models in the WC9 parameter domain, including here 
the extended C 11 model atom. Compared to the earlier subtypes, 
the WC9 stars have slower winds, as indicated by the smaller 
width of their lines. We therefore calculated one WC9 grid with 
iiaa = 1600 km/s and another with only !;„ = 1000 km/s. 

The two WO-type stars in our sample turned out to be very 
hot (T* x 200 kK). We therefore had to include higher ions of 
the iron-group elements, up to Fe xv, in their models. The ex- 
tremely broad emission lines indicate wind velocities as high as 
Voa x 5000 km/s, which is also much higher than for the WC 
stars. Hence, a special set of model calculations hat to be used to 
analyze the WO-type stars. 

The synthetic spectra of the main WC grid and the WC9 sub- 
grids will be made available from our websit43 simultaneously 
with the release of this paper. 

4. The analysis method 

The stars are analyzed by finding the PoWR model that most 
closely fits the observations. The first part of the fit procedure 
refers to the normalized line spectrum, and the second part to 
the fit of the spectral energy distribution (SED). Since WC spec- 
tra are full of lines, a rectification "by eye" can be ambiguous. 
Fortunately, most of the spectra we used are calibrated in abso- 
lute fluxes. We therefore rectify the spectra by dividing them by 
the model continuum. Consequently, the two steps (a) line fit and 
(b) SED fit are actually coupled and have to be iterated. 

The WC stars are divided into subtypes based on the ratios 
of the equivalent widths of certain emission lines. The stan- 

1 http : //www . astro . physik . uni -potsdam . de/PoWR . html 
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Fig. 1. Contours of constant equivalent widths of the WC emis- 
sion lines Cm 5696 A, Civ 5808 A and He 1 5876 A. The small 
dots represent models of the WC grid. Labels give absolute value 
of the equivalent width in A. 

dard c lassification system was summarized by Ivan der Huchtl 
(2001) and uses the ratio of Civ to Cm as a primary criterion. 
The main classification lines are Civ 5808 A and Cm 5696 A, 
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Fig. 2. Same as Fig.Q] but for the lines Ov 5590 A, Civ 5470 A 
and He n 5412 A. 



where Civ 5808 A is either one of two lines, together with the 
C m/C iv/He n-blend around 4650 A, that are extremely strong in 
the optical spectrum of early WC subtypes. 

One of our criteria for choosing the best-fit model is to en- 
sure in particular that we can fit these strongest lines. However, 
the enormous peak heights of in particular the earlier sub- 



Galactic WC stars 

Table 2. Model grid parameters 



WC grid WC9 grids 




55% 


Xc 


40% 


X °a 


5% 




0.16% 


log L/L Q 


5.3 


d„ /km s _I 


2000 1000/1600* 


D 


10 


Number of Levels 


He i 


17 


Hen 


16 


He m 


1 


Cn 


3 29 


Cm 


40 


Civ 


19 


Cv 


1 


On 


3 


Om 


33 


Oiv 


25 


Ov 


36 


Ovi 


15 


Ovn 


1 


Fe m-x" 


72 superlevels 



Notes. <a) Generic element, representing also Sc, Ti, V, Cr, Mn, Co, and 
Ni, with relative abundances to Fe as described in Graf ener et id] {20021 
{b) We calculated two WC9 grids with different velocities 

types can only be reached in the models at the price that other 
emission lines become stronger than observed. Thus, we of- 
ten had to make some compromise with respect to the over- 
all fit quality, and also paid attention to other significant lines 
such as O v 5590 A, the "diagnostic line pair" (Hen 5412 A and 
Civ 5470 A), and the carbon lines in the UV. For some stars in 
our sample, the red part of the optical spectrum was not avail- 
able, and our fit was restricted to the 4650A-blend and neigh- 
boring lines such as C iv 4441 A and sometimes O iv 3412 A. 

Hence, the line fit yields the stellar temperature T» and the 
"transformed radius" R t as main parameters, together with con- 
straints on the terminal wind velocity ««, and the chemical com- 
position. Fig. [3] shows as an example the observed optical spec- 
trum of the WC6 star WR 13, together with the synthetic spec- 
trum from the best-fit model. 

The PoWR models do not only provide a normalized spec- 
trum but also the absolute flux values of a star that can be used to 
fit the SED if flux-calibrated spectra and/or photometric values 
are available. The dilution of the model flux depends on the dis- 
tance and the reddening. For stars with a known distance modu- 
lus (D.M.), we can obtain the absolute magnitude M„ by 

M v = m v - D.M. - A v . (3) 

We note that the color indice s (v and b) ref er to the monochro- 
matic magnitudes defined bv lSrruthl (ll9681. The value of A v is 
provided by the SED fit where the reddening parameter Eb- V is a 
free parameter, as well as the detailed reddening law. Wheneve r 
adequate, we applied the reddening law from ISeatonl dl979h . 
which implies a fixed ratio of extinction Ay to color excess Eb-v 
of Ry = 3.1. If this turned o ut to be insufficient i n r eproducing 
the SED, we used the laws o flCardelli et all (1 19891) or lFitzpatrickl 
(Il999l) . which treat Ry as a free parameter, and adjusted the lat- 
ter to optimize the SED fit. (The reddening law and Ry are both 
indicated in the online material plots.) If mid-IR spectra were 
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Fig. 3. Optical spectrum of the WC6 star WR 13. The solid thin line is the observed spectrum, the thick dotted line represents the 
best-fit WC grid model. The primary model parameters are T* = 79.4 kK and log/? t //? = 0.5. The same model fits to most of 
the Galactic WC6 single stars, although it tends to underestimate the peak heights of Civ 5808 A and the Cra/Civ blend around 
4650 A. 



available for the S EP fit, we always chose the reddening law of 
iFitzpatrickl (1 19991) . 

An example SED fit is shown in Fig. [4] As the distance for 
this star is not known, an absolute magnitude M v = -4.42 mag 
was adopted from our subtype calibration. With the help of 
Fig.|4j we found that the luminosity of the grid model must be 
scaled to log L/L Q = 5.45 in order to fit the observed SED. 
Simultaneously, the color excess £V„ is also properly adjusted, 
together with the choice of the reddening law ("cardelli") and 
its parameter Ry- The fit implies a distance modulus of D.M. = 
12.15 mag. 

The color excess E/,- v in the Smith system is related to the 
more common Eb-v value from the Johnson system by Eb-v - 
1.21Ef,- v , while the extinction A„ is given by 

A v = (1.2W v + 036)E b - v (4) 

as described by Lundstrom & Stenholmj d!984l) . 

In general, it is hard to determine the distances of individual 
stars in our Galaxy, hence only for a subset of our sample are 
the distances known based on their membership to open clusters 
or associations (see Sect. [A}. For all other stars we follow the 
classical approach of assuming that all stars of the same subtype 
have the same absolute visual brightness M v . 

From the WC stars of our sample with known distances, the 
M„ values are shown in Fig. [5] The average relation is also in- 



Table 3. Calibration of absolute magnitudes for the different WC 
subtypes, as derived from stars with distances known from their 
cluster or association membership (cf. Fig.|5]l 



Subtype 


M v [mag] 


WC4 


-3.34 


WC5 


-4.12 


WC6 


-4.42 


WC7 


-4.18 


WC8 


-4.48 


WC9 


-5.13 


WN4/WCE 


-3.84" 


WN7/WCE 


-5.67° 



Notes. <n) Adopted WN calibration from lHamann et al.l d2006h 



dicated and tabulated in Table|3] It is unclear whether the as- 
sumption of a unique absolute brightness within each subtype 
is justified. The scatter in the M„ values within one subtype in 
Fig.|5]may only reflect the errors in the adopted stellar distances, 
but could also by partly intrinsic to the stars. Within these un- 
certainties, the resulting M v values increase towards later sub- 
types. This is similar to the c orresponding subtype calibration of 
the hydrogen-free WN stars (lHamann et al.ll2006h . We note that 
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log Ay A 

Fig. 4. Spectral energy distribution (SEP ) for the WC6 star WR5. Flux-calibrated observations (blue thin and noisy line) are from 
IUE in the UV. [Torres & Massevl d 19871) in the optical, and 2MASS (J, H, and K photometry) in the near-infrared. The synthetic 
spectrum is plotted as a red thick line. The continuum-only model flux is also shown for comparison. The model that fits the 
normalized line spectrum best has been selected from the WC standard grid (T* = 19 A kK, log Rt/R e = 0.5). 



even the same visual brightness would still imply a sharp depen- 
dence of the luminosity on subtype, because of their different 
bolometric corrections (see below). 
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Fig. 5. M v calibration for the different WC subtypes, using stars 
for which the distance is known from their membership to an 
open cluster or association (labels: WR number). The symbols 
connected by the dark line show the calculated average values 
for M v . The obtained values are listed in Table|3] 



For the two WN/WC stars without known distances, we 
adopt the M u c alibration for the corr esponding hydrogen-free 
WN subtype of Hamann et al. (2006), because their chemical 
composition is more similar to WN than WC stars. 



5. Results 

5.1. Stellar parameters 

The results of the analyses are compiled in Table|4]for our whole 
sample, now sorted by spectral subtype. The spectral fits and 
comments on individual stars are given in the appendix (online 
version only, see Sect. ?? andlAl. 

We begin by discussing the results from the line fit, chiefly 
the stellar temperature and the transformed radius. Figure[6] 
shows the locations of the analyzed stars in the log/? t -log T*- 
plane. The striking result visible in that diagram is the clear, 



one-dimensional sequence of the WC subtypes. They are aligned 
along the linear relation 

log/? t oc-21ogr t (5) 

with little scatter. The stellar temperature of the WC stars corre- 
lates tightly with their subtype. The two WO stars in our sample, 
however, do not appear to follow an extrapolation of the WC 
sequence, but are much hotter. 

Another group of points in Fig. [6] refers to stars that are 
known or highly suspected to be binaries with a luminous com- 
panion. We analyzed their spectra as if they were single stars. 
The parameters obtained by these "pseudo fits" clearly distin- 
guish the binaries from the sequence of single WC stars. This 
independently confirms the composite nature of their spectra. 
Thus, the "diluted emission line" (d.e.l.) criterion, which is in - 
sufficient for WN stars as emphasized by Hama nrTet alJ (|2006), 
seems to work fine for WC stars when distinguishing single stars 
from binaries. We conclude that those stars that follow the WC 
sequence in Fig.[6]do not suffer significantly from the line dilu- 
tion effect, i.e. they do not have a luminous companion. 

The tight correlation of the parameters with WC sub- 
type visible in Fig. [6] is now much clearer than in the first 
comprehensive analyses of Galactic WC stars by our group 
dKoesterke & Hamannl H995). Moreover, the subclass of weak- 
lined WC stars (WC-w) that we had introduced in the old paper 
is no longer necessary: all WC-w stars have since been identified 
as binaries. 

Some deviation from the straight alignment of subtypes is 
seen in Fig. [6] for the WC9 stars. The WC9 stars are character- 
ized by the appearance of strong C n-lines, and are significantly 
cooler than the other WC subtypes. 

In most of the WC9 spectra, the lines are narrower than in the 
simulations of our standard WC grid, which is calculated with a 
terminal wind speed of Voo = 2000 km s -1 . We therefore calcu- 
lated special grids in this parameter range with u OT = 1600 kms -1 
and = 1000km s _1 . 

The fit quality for the WC9 spectra differs from those of 
the other subtypes. While the "IUE long" UV spectra (1900- 
3300 A) are usually accurately reproduced and the "IUE short" 
range (1200-2000 A) is hardly met, the opposite is often true for 
the few WC9 stars with available UV spectra (WR69, WR92, 
WR 103). 
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In the optical part of the spectrum, most lines are repro- 
duced consistently by our models, although the peak heights 
of C in 4648 A and some lines between 6500 A and 7200 A are 
not recreated. The temperature steps in our grid are a bit too 
coarse to find the optimum fit in some cases, e.g. when the lines 
Civ 5808 A and He 1 5876 A are observed with almost the same 
peak heights. 

The WC9 stars of our sample cluster into two distinct groups 
in the log7? t -log 7Vplane, one cooler group around T* ^40kK 
and log R t = 1.0, which we refer to as group I, and another one 
with higher T* * 45 kK and a denser wind (log R t = 0.8), which 
we call group II. Could it be that group I consists of yet unde- 
tected binaries, while only the stars in group II are really single 
stars? 

In group I, only WR 104 (the "pinwheel") is an established 
binary, but all others are suspected to be binaries (see appendix 
Sect.|A]for corresponding details about WR59, WR65, WR69, 
and WR 88). On the other hand, the WC9d sta r WR 119 from 
the ho tter group II was checked for binarity by Williams et al. 
(2005) with a negative result. The same holds for WR 117, which 
is even hotter than both groups and closer to the WC8 stars in 
Fig lSl Two of the three WC9 stars in our sample that do not 
show dust (WR 81, and WR92) belong to the hotter group II, 
which seem to consist of the definitive non-binaries. The third 
one (WR88) seems to be a special case as discussed in Sect.lAl 

For very dense winds (log R t /Ro 55 0.4), the spectra no 
longer depend on the two indepe ndent paramete rs T„ and R t , 
but only on the product R t T^ (Hamann et al. 2003). The physical 
reason for this approximate parameter degeneracy is that in very 
dense winds, all radiation (including the continuum) emerges 
from layers that expand with nearly the terminal wind speed. 
Hence, under these conditions, the mass-loss rate is the only rel- 
evant parameter when comparing models with same luminos- 
ity and v m . Lines of constant R t T^ (or, equivalently, constant 
L/M 4 ^) are indicated by thin parallel lines in the lower part of 
Fig-El The model spectra in this range are nearly the same along 
these lines, hence the specific combination of R t and T% is not 
well-constrained. In other words, we see only rapidly moving 
layers of the atmosphere, while the slower part of the wind is 
opaque at all wavelengths, and therefore the radius of the hydro- 
static stellar core is not observable. This degeneracy primarily 
affects the WC4 stars. 

When available and not in conflict with our spectral fits, 
terminal ve l ocitie s co mpiled in TablelH were collecte d from 
iPrinia et al.l(ll990h and lNiedzielski & Sk"orzvnskild2002l) . These 
values had been derived from UV spectra, using the blue end of 
P Cygni line profiles. Therefore, they might differ from the ter- 
minal velocities obtained from optical emission lines. In these 
cases and for all remaining stars - mainly those without existing 
UV observations - Voo was roughly estimated by ourselves from 
the widths of the optical lines. We note that the models used for 
the line fits shown in the Appendix are not calculated with the in- 
dividual wind velocity of the star, but are taken from the models 
grids with fixed o„. The individual values are inserted, however, 
when Eq. is applied. 

We have so far presented those results derived from fits of 
the normalized line spectra. We now consider parameters which 
involve absolute dimensions, namely R„, M, and L, and those 
required to fit the SED in terms of absolute fluxes. 

The color exce ss Eh- „ , for w hich we use the monochromatic 
colors defined bv lSmithl d 19681) . of the best SED fit is included 
in Table|4] Further details, i.e. the applied reddening law and its 



Rv parameter if applicable, are printed on the fit plot shown for 
each star in the online Appendix. 

As explained above in Sect.lU the SED fit relies on either 
the adopted distance or absolute brightness, as indicated by the 
direction of the little arrows between the respective columns in 
Tableg] 

The basic free parameter of the SED fit was the stellar lumi- 
nosity. The models, which were all calculated for a fixed lumi- 
nosity of log L/Lq = 5.3, were scaled to match the observation. 
The parameters L and T„ were used to derive the stellar radius 
from Eq. ([TJ, and from the "transformed radius" R t we obtained 
the mass-loss rate via Eq. (f2]) under the assumption of D = 10. 
The final values of R*, M, and log L are included in Table[4] 

From the luminosities and the stellar temperature, we can 
construct the empirical Hertzsprung-Russell diagram (HRD, 
Fig.|7j. The stars with distances known from cluster or associ- 
ation memberships are represented by larger symbols, because 
their luminosities are more trustworthy than those relying on the 
subtype calibration. 

As one can see, the WC stars are located between the hy- 
drogen and helium zero-age main-sequence (ZAMS), except for 
some of the WC stars that fall on the cooler side of the hydro- 
gen ZAMS. The two WO stars have temperatures above those 
that would be expected from the He-ZAMS. The WC subtypes 
form a sequence in the HRD, from the WC9 stars, which are 
the coolest stars with the lowest luminosities, to the WC4 stars, 
which are hotter and more luminous. It is unclear whether the 
scatter is intrinsic to the stars, or caused by wrong distances. 
Taking the luminosities at face value, they range from about 
log L/L G = 5.0 for the faintest WC9 star to about log L/L = 5.8 
for the brightest WC6 star and the two WOs. Two of the WN/WC 
transition type stars (WR26 and WR 145) have exceptionally 
high luminosities, which are rather typical of WN stars. The 
HRD positions are discussed further in Sect. [6] as regards stel- 
lar evolution. 

All our program stars show strong mass-loss. The empiri- 
cal mass-loss rates from Table|4] are plotted versus their corre- 
sponding luminosities in Fig. [8] and follow a tight linear correla- 
tion in the double-logarithmic diagram, except for the WO stars. 
This is unsurprising if we remember that the WC stars basically 
follow the proportionality R t oc T* (cf. Eq.|5j- If all stars had 
the same u M , this would yield the mass loss-luminosity relation 
M oc I?I A This slope with a p ower of 3/4 was already sug- 
gested by Gra fener et alj d!998l) . who combined analyses from 
WC stars in the LMC and the Galaxy. The actual trend in Fig. [8] 
is slightly steeper, because v x also increases towards the earlier 
subtypes. The least square fit to our data gives 

log— ^— = (0.83 ±0.11)log— -9.10 + 0.57 (6) 

M Q /yr L Q 

Our empi rical data confirm ama zingly well the relation sug- 
gested by Nug is & Lamersl ([2000) - cf. Eq. (21) therein - based 
on the empirical data available at that time, when we assume a 
heavy element mass fraction in their formula of 45% (thick blue 
line in Figure[8]). The value of 45% results from our model as- 
sumptions of 40% carbon and 5% oxygen. 

Two more columns are included in Table[4] The current stel- 
lar mass M is deduced from the luminosity . For WC stars, the 
corresponding relation from Langer (1989) was applied, while 
for the transition types we adopted the relation for hydrogen-free 
WN stars from the same reference. The last column in Table|4] 
lists the values for the wind efficiency parameter rj, defined as 
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Fig. 6. The results of the WC analyses in the log T* -log R t -plane, which corresponds to the basic model parameters. The WC 
stars form a one-dimensional sequence with a slight offset for the WC9 stars. Stars appearing above the sequence turn out to be 
bi naries. The symbol sh apes indicate the different WC subtypes. The numbers next to the symbols are the WR numbers according 
to Ivan der Huchll ( 120011) . Binary systems analyzed as a pseudo fit have WR numbers in brackets. Stars with persistent dust emission 
are marked with a superscript "d". The thin grey lines in the lower part indicate the thick wind regime where stars can be shifted 
along these lines without significant changes in the synthetic spectrum. 



the ratio of wind momentum Mvca and photon momentum L/c 
per time: 



The "single scattering limit", at which every photon is scat- 
tered exactly once to accelerate the wind, corresponds to 77 = 1 . 
Hence, each photon has to be scattered on average 77 times 
to dri ve the wind with the empir ical parameters given in the 
table. iGrafener & Hamannl d2005l) constructed a hydrodynam- 
ically self-consistent model for the WC5 star WR 111 that 
achieved 77 = 2.54, thus showing that efficiency values of a few 
can be reached theoretically when multiple-scattering effects are 
properly taken into account. However, it seems questionable 
whether values of 77 S 5 can be explained. In our results many 
stars have wind efficiencies of 7/ > 5, or even 77 > 10, especially 
among the early WC subtypes. Interestingly, the two WO stars 
have only moderate 77 values. 

Possible solutions to this wind-driving problem include: (a) 
multiple scattering is even more efficient than hitherto thought, 
(b) there is an additional, yet unidentified wind-driving mech- 
anism, and (c) our mass-loss rates are drastic overestimates be- 



cause the clumping factor is much higher than the value of D-10 
that we adopted. 

An individual and detailed determination of the chemical 
composition of our program stars is beyond the scope of the 
present paper. For the WC stars, we checked that the compo- 
sition adopted for our model grid (with 40% carbon, cf. Tabled 
is roughly adequate. This can best be seen from the neighbor- 
ing "diagnostic line pair" He 11 5412 A and Civ 5470 A. These 
two lines form in the same zone of the wind and display a very 
similar dependence on T* and R t , while the ratio of their line 
strengths depends sensitively on the C/He abundance ratio. We 
find that a mass fraction ratio C:He of 40:55 is adequate for all 
WC stars in our sample, inclu ding the WC9 subtypes. As sug- 
gested bv lCrowther et al.l d2006l) . there is no indication that WC9 
stars are chemically less evolved than the earlier subtypes. 

Overall, models with an oxygen mass fraction of 5% seem to 
be appropriate for WC stars. However, the situation for the indi- 
vidual oxygen-line fit differs. The O iv-line at 341 1 A is usually 
nicely reproduced, while the fit quality of the O v line at 5590 A, 
which is used for classification, varies even within the subtypes. 
For the WC4 stars WR52 and WR 144, the line is significantly 
stronger. Thus, we used a model with an enhanced oxygen frac- 
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Fig. 7. Hertzsprung-Russell diagram of the 
Galactic WC stars analyzed in this work. Stars 
with known distances are represented by larger 
symbols. The sy mbols are labeled with the WR 
catalog number ( Ivan der Huchtl 1200 lh of the 
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tion of 15%. This raises the question of whether there might be 
an oxygen trend throughout the WC subtype sequence. There 
are indeed some stars of the subtypes WC5 and WC6 where we 
see that O v 5590 A is significantly stronger than predicted by 
the model, but there are also several examples where this line is 
perfectly reproduced with no more than the standard 5% oxygen 
mass fraction. We therefore conclude that there is no general in- 
crease in oxygen throughout the WC sequence, although there 
may possibly be an increate for the WC4 subtype. 

The WO stars, however, require models with even more oxy- 
gen, roughly 30% (by mass) instead of the 5% chosen for the 
WC grid. Interestingly the carbon abundance does not increase 
but stays at 40%. The remaining approximately 30% of the mass 
consists of helium, as the fraction of iron group elements is also 
the same as in the WC models, namely 0.16%. 

The spectra of the three WN/WC transition-type stars in our 
sample could not be fitted with models of WC-type composi- 
tion. Their spectral appearance resembles more those of the WN 
than the WC stars. To simulate these spec tra, we started from 
our models of hydrogen-free WN stars (see Ham ann et al.l2 006). 
These models contain mostly (98%) helium, plus 1.5% N and 
only 0.01 % C. While this nitrogen abundance also complies with 
the N features in our WN/WC spectra, the carbon abundance has 
to be increased in order to reproduce the C lines. 

We found that the adequate carbon abundances differ from 
star to star in this very special group. In particular, WR58 
(WN4/WCE) and WR 145 (WN7/WCE) contain a small amount 
of carbon, namely 0.1% and 0.5% by mass, respectively. In con- 
trast, WR 26 (WN7/WCE) even requires models with a carbon 



mass fraction as high as 20% to reproduce its spectrum. We note 
that we have omitted oxygen in our WN/WC special models to 
accelerate our calculations, which seems justified because the 
observed WN/WC spectra show hardly any oxygen features. 

The WC/WN star WR 126 is unique to our sample. The only 
notable emission lines are the C-He-blend around 4650 A, the 
Civ5808A-line, and three minor lines. Two of them can eas- 
ily be associated with Hen. The peak around 7 100 A, however, 
does not match any line usually seen in WC spectra. However, 
WN stars display a blend of N iv lines in this place. For this 
reason, this star is classified as WC5/WN. For our fit, we used 
a nitrogen-free WC star model with 20% (i.e. reduced) carbon 
mass fraction and 5% oxygen. 

Hence, it seems that for the transition-type stars neither the 
carbon abundance nor the position in the log ^ t -log 7Vplane (see 
Fig-Eb correlates with the spectral subtype. 

Summarizing the results for our whole sample, we found that 
all parameters show a close correlation with the spectral subtype 
(except for the transition types to WN). In Tabled we compile 
the mean parameters for each subtype. These results are dis- 
cussed further in Sect.|6]as regards the stellar evolution. 

5.2. Error estimations 

As explained in Sect. 15.11 the spectroscopic parameters (T», R t ) 
were determined by selecting the best-fitting synthetic line spec- 
trum from our model grids. Those grids are spaced by 0.05 in 
logT, and 0.1 dex in R t , which corresponds to our experience 
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Table 4. Parameters of the Galactic single WC stars 
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The arrows indicate whether we infer the absolute magnitude (M„) from a known distance modulus D.M., or whether the distance is calculated 
from M„ as obtained by the subtype magnitude calibration (cf. Table[3j. 

that this is roughly the accuracy with which the optimum line fit 
can be identified. 

For very dense winds (log R t /R 5 0.4), all emergent radia- 
tion originates from the fast-moving parts of the stellar wind, as 
discussed in Sect. 15.11 Hence, in this regime of parameter degen- 
eracy only the product R t T^ is constrained by the line fit. 

While the spectroscopic parameters are derived from the 
(normalized) line spectrum, the luminosity is obtained by fitting 



the SED. For most of our program stars, the photometrically cal- 
ibrated observations cover a wide spectral range from the UV 
to the IR, hence the color excess Eb- V can be determined to an 
accuracy of ±0.02 mag. The uncertainty in the distance is larger, 
which enters the derived luminosity oc d 2 and the mass-loss rate 
oc (P^ 2 . For many of our stars we must rely on the assumption that 
the u magnitude is constant per subtype (Sect.|4|. This whole cal- 
ibration, in turn, is based on those stars that presumably belong 
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Fig. 8. Empirical mass-loss rates versus luminosity for the an- 
alyzed Galactic WC stars. The subtypes are distinguished by 
different symbols as explained in the inlet. The numbers beside 
the symbols identify each star by its number in the WR catalog 
(Ivan der Huchtl 120011) . The red-dashed red line gives the least 
square fit to these re sults. The full, shaded fin es represent the re- 
lations proposed by [Nugis & Lamersl (2000) for WC stars with 
40% and 60% carbon, respectively. 



Table 5. Comparison of stellar parameters with previous analy- 
ses 
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2000) 



2002) 



bv lHillier& Miller! dl999l) . lDessart et al.l d2000t) studied another 
four Galactic WC stars, two of them, WR 90 (WC7 ) and WR 135 
(WC8), being single stars. ICrowther et al.l d2006l) analyzed the 
WC9d star WR 103. Their results are compared in Table[5]with 
those obtained in this work ("t.w."). The agreement among the 
stellar temperatures is very good. We tend to obtain slightly 
higher luminosities by 0. 1 dex. The mass-loss rates have a con- 
siderable scatter; we note that they were derived with the same 
assumption for th e clumping paramete r (D — 10) as in this work. 

Interestingly, IDessart et al.1 d2000l) also concluded that there 
is n o correlation between stellar temperature and subtype, as in 
our lKoesterke & Hamannl(fl9 95) paper, and in sharp contrast to 
our present result. The reason is obviously that their study was 
also confused by binaries that have diluted emission lines. 



to a certain open cluster or association. Apart from the possibil- 
ity that the membership assignment may be erroneous in indi- 
vidual cases, the distance to the clusters or associations are also 
continue to be debated between different authors, as particularly 
mentioned in appendix Sect.fA] To our impression, the distance 
moduli remain uncertain by about 0.75 mag, leading to an error 
margin of +0.3 dex for the luminosity. 

Systematic errors of our analyses are extremely difficult to 
quantify. Atomic data are one source of uncertainties. A more 
important question is the degree to the basic model assumptions, 
such as spherical symmetry and homogeneity, are adequate for 
real stars. A pl ethora of recent work deals with clumping in hot 
star winds (cf. lHamann et al.ll2008l) . Clumping can bias in par- 
ticular the empirical mass-loss rates. 

5.3. Comparison with previous WC analyses 

A l arger sample of Galact i c WC stars were analyzed in our group 
by iKoesterke & Hamannl (1 19951) . yet with un-blanketed PoWR 
model atmospheres. Their results did not reveal the clear corre- 
lation between the subtypes and the stellar parameters that we 
see now. The former analysis was partly confused by a couple 
of objects that have since been identified as binaries. In the case 
of WN stars, the improvement from un-blanketed to iron-line- 
blanketed models le d to about lOkK higher stellar temperatures 
( lHamann et al.l 12006b . Surprisingly, this is not the case for the 
WC stars. The newly obtained mass-loss rates are lower, mainly 
because of the microclumping correction (cf. Sect.01. 

A couple of Galactic WC stars were analyzed more recently 
with the help of the line-blanket ed CMFGEN model atmosphere 
code by iHillier & Miller! (Il998l) . WR 1 1 1 (WC5) was analyzed 



6. Stellar evolution 

To discuss the evolutionary origin of the WC stars, we now com- 
pare their empirical positions in the HRD with those of their po- 
tential progenitors, the WN stars. The results of this work are 
shown in Fig. [9] to gether with t hose of the Galactic WN stars 
studied by Haman n et al.l d2006l) . The WN stars form two dis- 
tinct groups, the very luminous "late" subtypes (WNL) with 
a significant fraction of hydrogen in their atmosphere, and the 
less-luminous "early" subtypes (WNE), most of which are com- 
pletely hydrogen-free. 

Like the WNE stars, the WC stars group at the left, hotter 
side of the ZAMS, while the hydrogen-containing WNL stars 
lie on the cooler side of the main sequence. The WO stars are 
even hotter than stars on the helium zero-age main sequence 
would be. We note that the luminosity range for the WC stars 
of early subtypes (WC4-7) and two WO stars is almost the same 
as for the hydrogen-free WNE subclass. The latest WC subtypes, 
WC8-9, are less luminous than any WN star. The large group of 
WNL stars, which are very luminous (log L/L Q > 6), have no 
counterparts to similarly luminous WC stars. 

These results provide an interesting test case for stellar 
evolution models. The s tellar evolution tracks calculated by 
iMevnet & Maeden d2003l) from the Geneva group are the stan- 
dard reference for massive-star evolution. In Figs. [10] and ITTl 
we compare the HRD of our WC sample with the tracks of 
Mevnet & Maeder ( 2003]) in the different versions with and 
without accounting for rotation, respectively. According to the 
chemical composition at the surface, the tracks are drawn in dif- 
ferent styles: hydrogen-rich phases are represented by thin black 
lines. When the hydrogen mass fraction of the atmosphere drops 
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the conclusion is that WC stars mainly arise from initial masses 
lower than 40 M . 

Our comparison of effective temperatures has some caveats, 
since these depend on the choice of the reference radius. The 
evolutionary tracks are plotted over the stellar effective tem- 
perature that refers to the hydrostatic stellar core. Our stellar 
temperature T* refers to the radius of Rosseland optical depth 
20 (cf. Sect. [3]). For the adopted velocity (and implied density) 
structure, this is close to the hydrostatic radiu s for most of the 
best-fit models. How ever, there are indications dlshii et alJll999fc 
iPetrovic et al.l2 006) that WR stars may have very extended, sub- 
photospheric layers, i.e. the hydrostatic core may actually be 
much smaller and hotter than derived from the observation of 
the stellar atmosphere. Moreover, for stars in the domain of the 
parameter degeneracy (i.e. WC4-5) discussed in Sect.[5] the hy- 
drostatic stellar radius can only be estimated from a significant 
inward extrapolation of the velocity field into optically thick lay- 
ers, which can lead to either an over- or under-estimation of the 
corresponding effective temperature. 

The empirical HRD of the WC stars shows a temperature gap 
between the helium main-sequence and the WC positions. While 
the WO subtypes are located to the left of the He-ZAMS (pre- 
sumably indicating the position of a hypothetical CO-ZAMS), 
the later WC subtypes are progressively cooler. The later the WC 
subtype, the lower the luminosity (implying a lower current mass 
according to the mass-luminosity relation), and the larger must 
be the sub-photospheric layer enhancing the stellar radius on top 
of the hydrostatic core. 



below 40% as in the WNL stars, the track is shown as a thick 
gray line. The hydrogen-free WN stage (Xu < 5%) is indicated 
by a thick dark line. In the WC stage (thick black/blue line), car- 
bon reaches more than 2 0% in the stellar enve lope. We chose this 
limit to be higher than in lHamann et al.l (120061) to ensure that the 
WC part of the tracks exclude any WN/WC transition stage. 

The evolutionary tracks predict that in the non-rotating sce- 
nario, massive stars with an initial mass of 60 M and above un- 
dergo a luminous blue variable (LBV) phase before entering the 
WNL stage. However, in the scenarios with rotation these stars 
skip the LBV p hase and enter the WNL stage directly after the 
main sequence (Mevnet & Maede rl2003h . Afterwards, the WNL 
stars evolve to hydrogen-free WNE stars at slightly lower lu- 
minosities, depending on the initial mass. They eventually enter 
the WC stage at roughly the helium main-sequence. According 
to the continuous mass-loss and the mass-luminosity relation, 
their luminosity decreases by a few tenths in log L before the 
gravitational core-collapse. 

Comparing now the Geneva tracks with our empirical HRD 
(Figs.fTOlandfTTli. we find relatively good agreement for the two 
WO stars of our sample. In particular the tracks without rota- 
tion and for initial masses > 60 M reach very hot and luminous 
endpoints prior to the supernova (SN) explosion. 

For the WC stars, however, the comparison ends less favor- 
ably. The first problem is posed by the luminosities. The evolu- 
tionary calculations without rotation (Fig.fTTIi predict that only 
tracks for initial masses above 37 M e return from the red super- 
giant (RSG) stage to the blue side of the HRD. The implied lumi- 
nosities are higher than those inferred from our WC sample. The 
tracks with rotation (Fig.flOb predict WR stars for initial masses 
above 22 M Q . Hence, the tracks are at least compatible with the 
most luminous WC stars of our sample, the WC6 subtypes. The 
bulk of WC stars, however, are less luminous than predicted by 
any of the tracks in their post-RSG phase. If we apply these re- 
lations between mass and luminosity from the Geneva tracks, 
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Fig. 1 0. HRD with the W C star positions and the tracks from 
iMevnet & Maederl d2003l) with an initial rotational velocity of 
u ro t = 300 km/s. The thick lines indicate the WR phases of the 
tracks. 

One could argue that all luminosity values are affected by 
(almost) all distance values for Galactic WR stars being inferred 
from cluster memberships. Even if some of the distances for in- 
dividual stars were wrong due to a false assignment to an asso- 
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Fig. 1 1. HRD with the W C star positions and the tracks from 
Mev net & Maederl d2003l) without rotation. The thick lines indi- 
cate the WR phases of the tracks. 



ciation or cluster, only their HRD positions would change while 
the overall picture would not be drastically affected, as long as 
cluster distances are not systematically underestimated. 

In contrast t o the Geneva stellar evolution models, 
Vanbever en et al.l JT998) showed that it is indeed possible to 
reach the WR stage for stars with a minimum initial mass of 
only 15 M Q , if the mass-loss rate during the r ed super-giant 
(RSG) stage is higher than given by the formula o flde Jager et al.l 
( 1988) adopted for the Geneva tracks. When adopting a higher 
RSG mass-loss, the resulting tracks (see Fig.fTZli predict, for in- 
stance, a luminosity of log L/L G « 5.2 for a star with Mj n ; t = 
20 M Q entering the post-RSG stage, and a final luminosity of 
logL/L =s 4.8 before explosion. Although in the models of 
IVanbeveren et al.l ([1998 ), WC phases only occur for stars with 
Mi n it > 25 Mq, the tracks cover almost the whole luminosity 
range of our WC sample. As for the results from the Geneva 
tracks, there is a discrepancy in the temperature T„ between the 
tracks and the evolutionary models. To date, none of the evolu- 
tionary calculations performed have been able to fully explain 
the obtained HRD positions, as they fail to reproduce the tem- 
peratures but provide enough indications that WC stars originate 
from a lower mass range. Both the Geneva and the Vanbeveren 
models imply that th e WC phase is close to the He-ZAMS. 
Graf ener et al.l (1201 2|) showed that extended sub-photospheric 
layers are not limited to hydrogen-rich stars on the upper part 
of the main sequence, but also occur for H-free stars on the He- 
ZAMS. This happens even in a lower mass regime if the clump- 
ing factor is sufficiently high. As briefly discussed above, such 
an envelope inflation would strongly reduce the effective tem- 
peratures of the stars, placing them exactly in the temperature 
regime we obtained in our work. Inflation might be a key factor 
in solving the temperature discrepancy. 

Another crucial aspect is the highly uncertain mass-los s 
rates during the RSG stage pointed out bv lMevnet et al.l (1201 ll) . 
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Fig. 12. H RD with the WC s tar p ositions and the tracks (thick 
lines) from lVanbeveren et al.l d!998l) using higher mass-loss rates 
during the RSG phase. The track for M m \ t = 25 M© is shown in 
light grey for clarity. In the WC phase, all tracks are drawn as 
thick dark/blue lines. 



which is partly caused by dust. For objects covered with dust, 
Ivan Loon et al.l (2005) discovered mass-loss rates that are sig- 
nificantly higher than those of visually bright ones. Such high 
rates would lead to stars evolving back from the red supergiant 
stage to the blue side of the HRD instead of exploding as a type 
II-P supernova. This would be in line with the observational 
constr aints for type II-P supernovae provided by ISmartt et al.l 
(2009), who found several SN II-P progenitors to be red su- 
pergiants, but none to have an initial mass of 18 M Q or higher. 
However, as there are known RSGs at least in the mass range 
between 18 M G and 25 M©, it is likely that they do not explode, 
but instead evolve to blue supergiants and eventually WR stars. 
Blue supergiants that have previously undergone an RSG stage 
would be He-enriche d. Objects with such an appearance have 
ind eed been found bylPrzv billa et al] (l2010h . 

IVanbeveren et al.l ([2007) also compared the parameters of 
the G alactic WN stars from our previous study (Haman n et al.l 
2006) with his tracks, which incorporate a higher RSG mass- 
loss rate (but no rotationally induced mixing), and concluded 
that their own tracks are in much better agreement with our re- 
sults than the Geneva tracks. 

Tracks that lead to WC stars with lower ini tial masses ques- 
tion the evolutio nary scenario s of IContil (119791) and also the re- 
vised version of lLanger etail (11994 that WNL stars evolve to 
WNE stars. Th e clear separat i on of WNL and WNE stars in the 
HRD found bv lHamann et all d2006l) suggests instead that WNL 
stars may not manage to get rid of their hydrogen layers, but ex- 
plode as type II supernovae either directly after their WNL phase 
or, more likely, aft er excursions into the LBV domain. More 
recen t results (e.g. iLiermann et alJl2010t ICrowther & Walbornl 
HoTD) indeed often discuss or assume that hydrogen-rich WNL 
(also labeled WNh) stars are core-hydrogen burning and cover 
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the high-mass end of the main sequence. Simply put, with in- 
creasing initial mass, the spectral types at the high-mass end of 
the main sequence change from Of to Of/WN and finally WNL. 

Moreover, our conclusions as regards stellar evolution can 
be explained in c onsistence with th e interpretation of super- 
novae statistics by Smit h et ail (1201 lh . They analyzed the differ- 
ent fractions of core-collapse supernovae (CCSN) using the data 
from the Lick Observatory Supernovae Search (LOSS), cover- 
ing mostly host galaxies with metallicities similar to Galactic 
metallicities (0.5-2Z o ). If the WNL stars do explode as super- 
novae without losing their hydrogen envelope, they will produce 
a type II supernova instead of the lb and Ic supernovae that are 
usually considered to be the end of WR evolution, because there 
is no hydrogen at a ll left in the WC stars and little in the WNE 
stars. ISntithT (12008) discussed how type Iln supernovae require 
stars with M m \ t > 50 M and these high mass stars can indeed 
explode without losing their hydrogen envelope. As type Iln su- 
pernovae require huge eruptive mass-loss that exceeds even the 
high WR mass -loss rates, the prima ry candidates are LB Vs, as 
suggested by iGal-Yam et al.l d2007l) a nd eventually proven for 
SN2005gl dGal-Yam & Leonardl2009h. The HRD posit ions ob- 
tained for the Galactic WNL stars dHamann et al]|2006l) corrob- 
orate our findings that the hydrogen-rich WNL stars are LBV 
progen itors which has also been suggested by ISmith & Contil 
(2008). All of the analyzed Galactic WNL stars seem to have 
initial masses higher than 60 M G . Their evolutionary sequence 
in the high mass range should therefore be 

Of -> Of /WNL <-> LBV WNL H - P oor] -» SN Iln. (8) 

There is an ongoing controversial debate about whether 
LBVs can be direct SN progenitors. In general, evolutionary cal- 
culations do not predict such a direct connection. The LBVs 
are characterized by their variability, which one cannot di- 
rectly measure after the star has explod e d, as e.g. in the case 
of SN2005gl dGal-Yam & Leonardll2009h . lKotak & Viiikl (120061) 
suggested LBVs as progenitors of certain supernovae due to sig- 
natures of progenitor variat ions in the SN spectrum. For SN 
2005gi. lTrundle et a l. (2008) argued that line profiles similar to 
the observed ones could only be found in LBVs. Finally, lumi- 
nosity and color, if these data exist, may indicate whether a SN 
progenitor was a blue star beyond the Humphrey-Davidson limit. 

Since we kno w WNL stars with v ery different hydro- 
gen abundances (|Hamann et alj 1200 6: Marti ns et al.l [2008; 
iLiermann et al.ll2.010h . it is likely that there are multiple excur- 
sions between a quiet WNL stage and LBV eruptions, as indi- 
cate d by the double- si de arrow in Eq. (|S). 

lDwarkad"asl d201 ll) questioned the scenario of LBVs being 
direct SN progenitors. He argued that the deduced mass-loss 
rates are usually based on the assumption of r -2 -density stratifi- 
cations, which might not be correct. He further argues that there 
should be at least a short WR phase between an LBV outburst 
and a SN explosion. For a type Iln supernova, this phase would 
correspond to a WNL type with a low hydrogen abundance. 
We therefore add the suffix "H-poor" to this possible post-LBV 
W NL stage in Eq. d8]l 

Smith et al] d201 ll) discussed that there seem to be two ways 
for LBVs, either by exploding or turning into another WR stage 
to eventually explode as SN lb or Ic. Since the results of our WC 
analyses suggest that WNE and WC stars have progenitors of 
lower mass, this decision conflict would be solved. For an initial 
mass range from approximately 20 M to 50 M , the WR stars 
would arise as post-RSG stars with most of their hydrogen being 
lost during the RSG stage, possibly skipping a hydrogen-rich 



WN stage that is not observed for this mass range. The evolu- 
tionary sequence leading to WC stars might therefore be 

O RSG WNE -> WC -> SN Ib/c. (9) 

However, there is an upper mass limit for red supergiants. 
Although the exact value is un clear, it might be somewhere be- 
tween 25 M Q and 50 M Q (e.g. iMevnet & Maederi 120031. ISmarttl 
2009), meaning that stars with higher initial masses do not be- 
come RSGs. Stars in the range above the RSG mass limit but 
below the WNL mass range might therefore undergo an LBV 
phase, but do not appear as WNL stars at any earlier stage. In ad- 
dition, these "lower mass LBVs" do not seem to explode as type 
II supernovae but instead lose their hydrogen layers and eventu- 
ally become hydrogen-poor stars, such as WNE, WC, and WO 
stars. Assuming a mass limit for red supergiants significantly 
below 50 M , the evolutionary sequence for the WC stars above 
would be 

O LBV -» WNE — > WC -> SN Ib/c. (10) 

Our empirical HRD positions for the Galactic WC stars in- 
stead support the post-RSG origin of Eq. (O. All WC stars are 
located below the 40 M tracks, independently of their detailed 
scenario assumptions. Even if the RSG upper mass limit is some- 
what lower, most if not all WC stars of our sample must be at- 
tributed to the post-RSG evolution. 

Regardless of how the WC stars are produced, if they ex- 
plode as a SN, it clearly has to be of type lb or Ic, as these stars 
are hydrogen-free. While there are known progenitor stars for 
type II supernovae, there are none for type lb and Ic. Therefore 
it is not yet clear if WNE or WC stars will undergo a supernova 
explosion at all. Instead they might just collapse to a black hole 
without a high luminosity outburst. 

To support the mass ranges for t his new sequences , we take 
a closer look at the discussions from lSmith et al.l d201 lh . assum- 
ing that WC stars always end up in a supernova. They tried 
to assign the different core collapse supernova types to mass 
ranges, so that the fractions fit with a Salpeter IMF. With the 
assumption all supernovae arise from single-star evolution they 
end up with a maximum mass of 22 M Q for type II supernovae, 
even including the SNelln here. That would imply that all stars 
above 22 M would completely lose their hydrogen, probably 
becoming a WNE and finally a WC or WO star. Although this 
range would fit with the lower end of our WC initial mass range, 
this scenario ignores the previously suggested WNL evolution 
and contradicts not only the idea, but also the observation of a 
SNe Iln with a high-mass progenitor. 

Even though they are only referring to core-collapse super- 
novae, it is clear that binaries have to be taken into account. In 
contrast to the assumption that only single-star evolution domi- 
nates the SNe statistic, one could assume that binary evolution 
is the only possible way for type Ib/c SNe. That would of course 
assume that all WNE, WC, and WO stars are binary systems or 
have at least been formed by close binary evolution. Although a 
number of systems might have formed in this way, it is highly 
unlikely that all WR stars (except WNL stars) share this origin. 
Stellar evolution models might not cover the right mass range 
for these stars yet, but they do show that WNE and WC stars can 
be formed without req uiring close binary evolution. 

Smit h et alj d201 lh therefore presented two ad hoc scenarios 
combining the ideas of single star and binary evolution. A frac- 
tion of 8.8% of their SN sample are of type Iln while a total 
fraction of 26% are type I CCSNe (hereafter labeled SNe Ib/c). 
Among these 26%, some originate from single stars and others 
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from binaries. As an ad hoc approach, they assumed that the 
single star fraction should be as large as the SNe Iln fraction, 
namely 8.8%. The remaining fraction is considered as a result 
of close binary evolution, together with the SNe lib, which are 
actually type lb SNe with small amounts of hydrogen visible 
only in the early stage of the SN. For the type II-P and II-L SNe 
associated with the low initial masses, they ended up with an 
upper limit to M; n i t of 23.1 M Q , Above this initial mass, all sin- 
gle stars explode as either type Ib/c or Iln supernovae. Were a 
strict mass limit assumed to separate the SNe Ib/c and Iln sce- 
narios, thjs_wou!dbeMini t = 37 M under the assumptions stated 
above. [Smith et al.1 d201 lb assigned the SNe Iln to the lower mass 
range between 23 M and 37 M , but according to the results of 
this work and the observation of iGal-Yam & Leon ard (2009|), it 
seems likely that the converse is true. The Iln might then belong 
to stars with M- m \ t > 37 M , possibly following the evolutionary 
scenario ([8]), and the SNe Ib/c come from WC stars with initial 
masses between 23 and 37 M , which is not so far from our pro- 
posal based on the HRD positions relative to the evolutionary 
tracks. 

The limits of the mass ranges are of course only rough es- 
timations. Neither our track comparisons, nor the ad hoc as- 
sumption of the single star fraction in the SNe Ib/c fraction are 
based on detailed calculations. The limits therefore can change 
in both directions, e.g. depending on the values of the RSG mass- 
loss rates and the number of binaries causing type Ib/c SNe. A 
larger single-star fraction in the ad hoc assumption would de- 
crease the minimum initial mass, while a lower RSG mass loss 
woul d raise it. A s an ad ditional channel leading to SNe of type 
Ib/c, ISmith et all d201ll) suggested underluminous H-poor stars. 
We consider this idea as speculative, because these objects have 
never been observed. Admittedly, if they formed by Roche-lobe 
overflow in close-binary systems, such stars might escape detec- 
tion. However, on the basis of the SN statistics, the latter case is 
covered by the binary fraction that was already subtracted from 
the number of Ib/c supernovae. Despite all the remaining un- 
certainties, the scenario demonstrates that the observed CCSN 
statistics in general can be explained in-line with our empirical 
results. 

We have so far neglected two stars that do not seem to follow 
one of the newly introduced scenarios ([8]l or (0, namely the two 
W02 stars WR 102 and WR 142. In contrast to the early WC 
subtypes they appear as hot, but also very lu minous stars. Their 
positio ns are close to the end of the tracks of iMevnet & Maederl 
(2003) with favoring the non-rotational models. Their positions 
therefore indicate that these two WO stars at least are stars with 
Afinit > 40 M Q . The end of each Geneva track is marked by a SN 
symbol, which is not entirely appropriate as the track calcula- 
tions end when there is no more helium left in the stellar core. It 
has been suggested that WO stars might alrea dy be in a state of 
carbon burning (e.g. Bar low & Hummerll982l) . The low number 
of known WO stars - just four of them in the Milky Way - and 
the high temperature of the corresponding PoWR model com- 
pared to the Geneva track endings support this idea. However, 
owing to the absence of unambiguously identified helium lines in 
the optical spectrum, our models might overestimate the helium 
abundance in the WO stars. Only if these stars were completely 
helium-free, we would have a clear evidence of a carbon-burning 
star. 

With the sequences ([8]l and (0, we have evolutionary sce- 
narios for the WC stars with masses of up to 40 M or maybe 
50 M and another for the WNL stars with M init > 60 M , sug- 
gested by the WNL positions and the SNIIn constraints. The 
observed WC luminosities also roughly agree with the tracks of 



IVanbeveren etal] {1998) for this mass range. The tracks do not 
reproduce the obtained temperatures but this might be due to 
envelope inflation, which has never been included in any evo- 
lutionary track to date. The WO stars might actually originate 
from stars with initial masses between the two mass bins. To 
end up as WO star, a star might fellow a scenario close to the 
one of lLangeretal] (11994 . which is also a kind of mixture of 
scenarios (0 and ([8} 

O -> WNL -> LBV -> WNE? -> WO -> SN Ib/c. (11) 

In contrast to the original assumption o flLanger et alJ d 19941) . the 
WC phase is skipped here. If the WO stars were to evolve from 
WC stars, we would expect to see WC stars with even higher lu- 
minosities than the WO stars, which are not observed. The situa- 
tion is not much better for WNE stars, hence this stage is tagged 
with a question mark in (fTTT i as there are no stars observed in 
a luminosity range that would correspond to the luminosities of 
the WO stars. However, this might be just a consequence of the 
low number statistics, as the WNE lifetimes are a factor of be- 
tween three and ten short er than those of the WC s tars, accord- 
ing to the calculations of IMevnet & Maederl d2003l) . In the case 
of He-burning WO stars, their lifetimes would be comparable to 
those of WC stars and with only two known W02 stars in the 
Milky Way, the absence of such WNE stars in the corresponding 
luminosity range would be merely a statistical effect. 

The final SN explosion, if there is one, should be of either 
type lb or Ic. The latter would assume that such stars have lost 
all their helium. This is not the case in our current WO models, 
which contain a helium mass fraction of 30%. However, as men- 
tioned above, we cannot rule out that these stars might already 
be helium-depleted. 

Alternatively, one could imagine that WO stars form from 
rapidly rotating WNL stars, with a very short or even without a 
WNE phase. This would be in line with our findings that both 
WO stars in our sample require a very high rotational velocity 
of 1000 km/s to reproduce their round-shaped spectral lines. On 
the basis of these assumptions, the evolutionary scenario would 
be 

O -> WNL -> WO -> GRB. (12) 

Stellar evolution calculations currently support neither sce- 
nario ([TTj) nor JT2l . The most recent Geneva models 
dMaeder & Mevnel201 ll) predict a spin-down during the main- 
sequence stage, leading to lower rotational velocities at higher 
masses. However, we do see extremely broad lines in the WO 
spectra that can only explained by very rapid rotation. We spec- 
ulate that these stars were either spun up by mass transfer from 
a companion, or that the transfer and loss of angular momen- 
tum remains inadequately described by the current evolutionary 
models. 

In total, the modified scenario of ISmith et alJ d201 ll) corrob- 
orates our conclusion that WC stars arise from those of lower 
initial masses than previously expected, as raised by the HRD 
positions. It is also in line with the assumption that the very 
luminous WNL stars are core hydrogen-burning stars that will 
never reach the WNE and WC stages. However, the SN scenario 
does not include metallicity effects and their binary assumptions 
are just ad-hoc. As the formation of WR stars might strongly de- 
pend on the metallicity of the host galaxy, the SN statistics show 
a plausible and interesting scenario that is in line with our results, 
but one that should not be taken as an independent argument for 
the WC mass range. 

Metallicity effects may be illustrated by comparing WC stars 
in the Galaxy and the Large Magellanic Cloud (LMC). The 
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six W C4 s tars in the LMC were analyzed by IGrafener et alJ 
(Il998h and ICrowther etaD d2002h . According to the latter pa- 
per, the_JurninosMes_ofJhese stars range from logL = 5.4 to 
5.7. IGrafener et al] dl998l) obtained slightly lower values (typi- 
cally by 0.2 dex), but used models without clumping and iron- 
line blanketing at that time. Compared to the average luminosity 
of the WC4 stars in our Galactic sample (see Tabled, it seems 
that their counterparts in the LMC are significantly brighter. The 
mass-loss rates, on the other hand, are about the same (the as- 
sumptions on clumping are identical in both works), which is 
obviously the reason why their line spectra are similar. One can 
conclude that (i) WC4 stars in the LMC have a higher current 
mass than Galactic stars of the same type and (ii) a higher lumi- 
nosity is needed at LMC metallicity to drive the same mass loss. 
For probably related reasons, WC stars of later subtypes and thus 
larger pseudo-photospheric radii do not exist in the LMC. 

7. Conclusions 

We have analyzed 56 Galactic WC stars by comparing their 
line spectra to simulations from a grid of line-blanketed PoWR 
model atmospheres to derive the stellar parameters. In addition, 
the SED was fitted by adjusting the reddening and scaling the lu- 
minosity to match flux-calibrated spectra and photometric mea- 
surements. For the first time, a clear correlation between the 
spectroscopic subtype and the stellar parameters (Tabled has 
been found. 

Our results also show that for single WC stars the subtypes 
are aligned along a strip in the log 7Vlog R t -plane with the WC9 
stars having a slight offset. Moreover, in the M-L-plane, we ob- 
tained an M-L-relation for WC stars in the form of M oc L 8 , 
which is in good agreement with Nugis & LamerJ (l2000h . 

For the chemical composition, the helium-to-carbon ratio 
seems to be rather uniform within the WC subclass (about 55:40 
by mass). The oxygen abundance is of the order of 5%, but might 
be slightly higher for the earliest subtypes and is about 30% for 
the WO stars. 

For the termi nal wind velocity, we used values from 
iPrinia et all d 19901) and iNiedzielski & Skorzvnskl d2002l) when 
available and scaled them for the remaining stars before calcu- 
lating the resulting stellar parameters. Table[6] reveals a signifi- 
cant increase in the mean value of v m per subtype. Starting with 
values of about v m w 1400 km/s for the WC9 stars with some 
of them having even smaller values (close to 1000 km/s), v m 
quickly increases for earlier subtypes, ending up as w 3300 km/s 
for WC4 s tars. The mass es obtained from the luminosity via the 
relation of lLangerl |l989) seem to roughly increase from late to 
early subtypes, but these values have to be interpreted with care. 
The underlying mass-luminosity relation is indeed independent 
of the scenario that led to the WR stage. 

The stars marked as WC9d display a significant excess even 
in the near-infrared, a finding that remains unexplained. Some of 
these stars might be part of a binary system. We obtained similar 
parameters as in our a binary pseudo fit of the well-known binary 
system W R 104 for some s tars th at are thought to be single stars. 
However, Williams et al.1 d2005h found indications that not all 
WC9d-stars might have a companion. 

The WO stars WR 102 and WR 142 (both of W02 subtype) 
differ significantly from the WC stars. They contain much more 
oxygen {Xo a 30%) and have both much higher temperatures 
(T t » 200 kK), and extremely fast winds with v m ~ 5000 km/s. 
They are very compact objects with masses higher than those 
of the WC stars and stellar radii smaller than 1 R & . Their round 
emission-line profiles indicate that these stars are rapid rotators. 



Table 6. Mean WC star parameters per subtype 



Subtype 


T, 
[kK] 


[km/s] 


R, 
[Re] 


log M" 
[M /yr] 


log L 
[Lei 


Mwc* 
[M Q ] 


WU2 


2U0 


5000 


0.0 


— j.Uo 


J. / 




WC4 


117 


3310 


1.0 


-4.65 


5.2 


10 


WC5 


83 


2780 


2.2 


-4.62 


5.3 


12 


WC6 


78 


2270 


2.9 


-4.61 


5.5 


14 


WC7 


71 


2010 


2.9 


-4.79 


5.3 


11 


WC8 


60 


1810 


4.2 


-4.80 


5.3 


11 


WC9 


44 


1390 


6.6 


-4.80 


5.2 


10 



Notes. <c,) Mass-loss rates are calculated with a volume-filling factor 
of / v = 0. 1 WC mass es are calculated from luminosities using the 
relation of Langed d!989l) 



The WN/WC stars are WN stars with an enhanced carbon 
fraction that have WN-type spectra. They are fitted with special 
models that reflect these compositions, revealing that their ba- 
sic parameters {T*, R ( ) are already close to those of the WC se- 
quence. The star WR 126 is sometimes thrown into this scheme 
because it is classified as WC/WN. In all cases, WR 126 does not 
fit into any of the former groups, and seems to be an exception 
that has either its own unique history or is a so far undetected 
binary system, this last possibility being supported by its prox- 
imity to the binary-pseudo-fit region in the log 7Vlog /? r plane. 

The obtained positions of the WC stars in the HRD chal- 
lenge the standard scenario that WC stars generally arise from 
very massive stars. Our results instead indicate that WC stars 
mostly originate from an intermediate-high mass range with ini- 
tial masses of between 20 M e and 40 M or maybe 50 M Q . This 
suggests that they are descendants of post-RSG or - for the 
higher masses - LBV stars, which have lost enough of their hy- 
drogen envelope that they directly enter the WNE stage before 
they eventually become WC stars. 

Although rotation decreases the minimum mass that is re- 
quired to reach the WR stage, it seems that the scenarios with- 
out initial rotation provide a closer fit to the obser vations. For 
the WC stars, the tracks of Vanb everen et al.l ([1998) did not in- 
clude rotation, but instead higher RSG mass-loss rates. These 
tracks are in much better agreement with our WC results than 
the Geneva tracks. The positions of the WO s tars in the HRD are 
closer to the non-rotational Geneva tracks of lMevneT & Maeder 
(120031 which appears to contradict the afore mentioned results 
for their line profiles, but it could simply indicate that fast rota- 
tion occurs only during the later stages of evolution. 

The evolutionary status of the WC9 stars is quite obscure, 
at least if they are single stars. They might originate from stars 
with initial masses of around 20 M Q , but their evolutionary fate is 
completely unclear. As a first approach, one could imagine that 
they evolve to earlier WC subtypes, but that would have to be 
consistent with their increasing luminosity, as for early WC sub- 
types the luminosities are much higher than for the WC9 stars. 
Such an increase in luminosity would contradict the current as- 
sumptions of the mass-luminosity relation that correlates larger 
masses with larger luminosities. 

The WNL stars, using the term for WN stars containing hy- 
drogen, seem to be different from the rest of the WR stars. Their 
HRD positions suggest that these stars might not actually be in 
the phase of helium burning traditionally expected of all WR 
stars. They might instead simply be core-hydrogen burning stars 
with very high masses and mass-loss rates, making them a kind 
of extension of the O stars at the high-mass end of the main se- 
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Table 7. Suggested single-star evolution scenarios based on WN 
and WC analyses 



Mm [M ] 


stellar evolution 


8-15 


OB^RSG^SNII-P 


15-20 


OB RSG BSG SNII-L 


20-45 


O RSG WNE ^ WC ^ SNIb/c 


45-60 


O WNL LBV/WNE? -> WO -> SNIb/c 


>60 


O -> Of/WNL <-> LBV [-> WNL] -> SNIIn 



quence. They are therefore sometimes not considered as "real" 
Wolf-Rayet stars. While this might seem logical in the way of 
using the term Wolf-Rayet stars for helium-burning stars, it is 
incorrect, as WR stars are not defined in terms of their core burn- 
ing status but only their spectral characteristics. Nevertheless, 
the differences in the evolutionary status between the WNL and 
the other WR stars should be carefully taken into account for 
future stellar evolution calculations. 

The WO stars, at least those of W02 subtype, seem to differ 
significantly from WC stars, as we conclude from the different 
parameters and HRD positions. However, their evolutionary sta- 
tus remains unclear. As a tentative explanation, we suggest the 
scenarios (fTTI) and lfl2b . With only four known WO stars in our 
Galaxy, the observational basis is much worse than for either 
WN or WC stars, especially as our results indicate that WO stars 
are not descendants of WC stars. Scenario (fTTT i probably leads 
only in special situations to a WO star, making WO stars such a 
rare subclass. 

When comparing the obtained HRD positions with evolu- 
tionary models, we see a significant temperature discrepancy for 
WC5 stars and all later subtypes. A promising candidate to solve 
this problem is enve lope inflation, as described for H-free stars in 
iGrafener etail (120121) . Unlike the later subtypes, the WC4 stars 
and also the very different WO stars do not seem to be inflated. 

A summary of the suggested single-star evolutionary scenar- 
ios for massive stars based on the discussion in Sect. [6] can be 
found in TableQ As we found the upper mass limit for the WC 
stars to be somewhere in-between 40 M e and 50 M , we simply 
give the mean value of 45 M , All of the mass limits specified 
in Table|7]are just rough suggestions, which are based on our re- 
sults and the referenced comparisons to SN statistics and stellar 
evolution models. 

Detailed stellar evolution models, including higher RSG 
mass-loss rates and rotation, would enable us to check whether 
our suggested evolutionary scenarios can be numerically repro- 
duced. More reliable distance determinations for the Galactic 
WR stars would help to verify the low luminosities obtained 
in this work. Detailed analyses of WC and WO stars in other 
galaxies would help us to circumvent the distance problems and 
enhance our proposed scenario in order to check how it depends 
on metallicity. 

Both the WC and WO star parameters can provide con- 
straints on the wind-driving mechanisms and the WC mass- 
loss rates can influence the chemical evolution of our Galaxy. 
Furthermore, WC and WO stars are keystones in the present the- 
ories for the evolution of massive stars. In the past years, our 
picture of massive stars has changed and new evolutionary sce- 
narios have been proposed. We hope that our results and sug- 
gestions can contribute to these discussions and help us to better 
understand the origin and evolution not only of the WC stars, but 
of Wolf-Rayet stars in general. 



PoWR model remark: The model grid used for the Galactic 
WC stars will be published on our PoWR website^] with the 
model files being ready-to-use. The interface will allow you to 
select and view the basic model parameters before downloading 
the model itself. 
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Appendix A: Comments on individual stars 

In this appendix we compile some additional information about 
individual stars of our sample, especially their distance and bi- 
nary status: 

W R4 is a WC5 star that is listed as SB1 in Ivan der Huchtl 
(200lJ) owing to its short-periodic photomet ric variability mea- 
sured bv lRustamov & Ch erepashchuk (1989), but it displays nei- 
ther dil uted emission nor a bsorption lines from a possible com- 
panion. [SmithetaLl d 19901) discussed whether WR4 belongs to 
an H i bubble, for which they estimate a distance of between 
1.6kpc and 2.73 kpc depending on the method. The smaller 
value is close to our spectroscopic distance of 1.64 kpc, based 
on the M„ calibration. 

WR8 is classified as WN7/WCE. Its spectrum resembles 
a WN7 t ype st ar, but with stronger carbon lines than usual. 
iNiemelal (Il99ll) measured large amplitude radial-velocity vari- 
ations with a possible period of 38.4 d, where the carbon and ni- 
trogen lines are in antiphase. This is indicative of a WN+WC 
binary system. We attempted to reproduce the spectrum us- 
ing a single-star model, which resembles a hydrogen-free WN 
star with e nhanced carbon, but failed to o btain a convincing fit 
(Fig. B.3). ILundstrom & Stenholml (1 19841) estimated the modu- 
lus to be D.M. = 12.7 mag based on a membership of the Anon 
Pup a association. 

WR14 is a WC7 star, for which photometric va riability 
has been reported by various authors, and discussed bv lShvlaia 
(1990) in terms of a possibl e compact companion . It displays 
non-thermal radio emission (Chapman et alT IT999l) . which can 
normally be attributed to colliding winds in binary systems. 
We found that its spectrum can be reproduced with a typi- 
cal WC model for single stars, and conclude that there is at 
least no bright companion th at dilutes the WC emission lines. 
ILundstrom & Stenholml d 19841) identified the star as a member 
of Anon Vel a, estimating a distance of D.M. = 11.5 mag. 

WR15 has the spectral type WC6. From the width of its 
lines, we estimated a terminal wind velocity of 2600 km s~', 
slightly higher values than typical for this subtype. The dis- 
tance of D.M. = 12.0 mag i s based on the possible membe rship 
of Anon Velb suggested by Lundstrom & S tenholml dl984l) . We 
prefer to used this distance method ove r that used to calculat e 
the spectroscopic parallaxes obtained bv lConti & Vaccal (Q990), 
who derived a value of D.M. = 1 1.28 mag for WR 15. 

WR23 is another a W C6 star and belongs to the Car 
OBI association, for which ISmithl ([2006) obtained a distance 
modulus of D.M. = 11.8 mag for the Homunculus nebula 
around 77 Car. Earlier distance calculations for Car OBI were 
P.M. = 12.55 mag dMassev & Johnson! 119931) and 12.1 mag 
dLundstrom & Stenholmll 19841) . 

WR26 is classified as WN7/WCE. The spectrum contains 
a very strong Civ 5808 A line, which requires a model with 
about a 20% carbon mass fraction to be reproduced. The fit to 
our single-star model is satisfactory. The high carbon abundance 
might indicate that WR 26 is undergoing a transition from a WN 
to a WC star. 

WR33 has the spectral type WC5. The UV spectrum is in- 
dicative of a high wind velocity of about 3000 km s" 1 , but other- 
wise looks normal. The optical spectrum (from lTorres & Massevl 
[T987h shows very broad emission features between 4000 A and 
4600 A that cannot be reproduced by any of our models. 

WR38 displays a WC4 spectrum. Our fit is not very 
consistent, possibly indicating that the lines are diluted by a 
companion's contin uum. From a weak cluster around WR 38, 
Shorl in et al.l d2004l) derived a very large distance of D.M. = 



15.8 mag, which would imply a very high luminosity of this star. 
In contrast, Wallace et al . (2005) obtained a distance modulus of 
15.0+j q mag from the HST photometry of the same cluster. Our 
photometric distance, based on the subtype calibration, yields 
D.M. = 15.1 mag, in perfect agreement with the latter work. 

WR39, classified as WC7+OB? dvan der Huchdl200lt) . is an 
interesting test case for the "diluted emission-line" binarity cri- 
terion. When we analyze its line spectrum assuming that the 
star is single, the obtained parameters place the star distinctly 
above the sequence of single WC stars in the log T*-logR t - 
plane (see Fig.[6]l. Moreover, W R39 shows non-thermal radio 
emission (Chapma n et al.ll 19991) . which is normally attributed to 
colliding winds in binary systems. We therefore conclude that 
WR 39 is indeed most likely a binary. 

WR45 is a WC6 star. Unfortunately, the only spectrum 
available to us is limited to 34 10-4730 A and does not cover 
many diagnostic lines. Therefore, the parameters we obtained 
for this star are uncertain. 

WR50 is a doub le-lined spectroscopic binary WC7+OB 
dvan d er I lucht 1200 If) with a light-curve period of Id 
dvan Genderen et all 1 19911) . We analyzed its spectrum as if it 
were from a single star, and obtained parameters that are obvi- 
ously affected by the d.e.l. effect (cf. Fig.[6]> . The distance mod- 
ulus is D.M. = 12.8 mag according to I Vazquez et al.l (|2005). 

WR52 is one of only five known Galactic WC4 stars, and 
does not show indications of binarity. The spectral fit Fig. B.16) 
requires models with higher oxygen abundances (10-15% by 
mass) than we used for our standard WC grid (5%). This agrees 
with our findings of highe r oxygen abundances for WC4 and 
WO stars. As discussed by Chu & Treffersl dl98ll) . a kinematic 
distance of 2 kpc and a photometric distance of 4 kpc were de- 
rived from the associated nebula. From our subtype calibration, 
we obtained a distance modulus of D.M. = 12.67 mag, corre- 
sponding to 3.4 kpc. 

WR53 is classified as WC8 with the "d" indicating persis- 
tent dust. In contrast to the WC9 stars, dust formation is uncom- 
mon for this subtype, and might indicate colliding winds. The fit 
to our models, calculated with the full C 11 i on, remains remark- 
ably poor (Fig. B.17) for unknown reasons. Ma rtin et alj (2007) 
found an expanding H 1 shell probably associated with this star, 
encircling an optical emission nebula, and estimated the distance 
of the shell to be 4+1 kpc. From our subtype calibration, we ob- 
tained 2.9 kpc. 

WR58 is a transition-type star, classified as WN4/WC. For 
this star, a carbon mass fraction of 0.1% is enough to reproduce 
the Civ 5808 A-line. 

WR5 9 is classified as WC9d. The star might be a binary, 
as IWilliams et alj (|2005) found Balmer absorption features by 
comparing the spectral lines to those of a non-dusty WC9 star. 

WR64 is classified as WC7 . The only spectrum available 
to us is that of lTorr es & Massev ( 1987) and covers only a short 
wavelength range, making our analysis less precise. 

WR65, cl assified as WC9d, shows vari able X-ray emission. 
According to Oskinova & Hamann (2008), the emission orig- 
inates from the wind-wind collision in a massive binary sys- 
tem, and the variability is caused by t he dif ferent absorption 
columns along the orbit. Williams et alJ d2005l) detected absorp- 
tion features from the Balmer lines of hydrogen. We analyzed 
the spectrum as if it were from a single WC star and found 
no indication of a composite nature. However, if W R 65 were a 
memb er of the CirOBl association as suggested bv lLortet et al.l 
dl987l) . it would be by far the brightest WC star in o ur sam- 
ple (M„ = -7.02 mag adopting a distance modulus from lTurnerl 
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(119961) of 12.57 mag). As this seems unlikely, we conclude that 
WR65 is located in the foreground of CirOBl, or it has a com- 
panion that contributes significantly to its total brightness. In any 
case, we cannot employ this star for our M„ versus subtype cali- 
bration. 

WR68 is of subty pe WC7 a nd a possible member of the 
CirOBl association. lTurneri dl996) deduced a distance of D.M. 
= 12.57 mag, based on the cluster Pismis 20, which is part of 
CirOBl. Adopting this number leads to a plausible luminos- 
ity and HRD position. However, the distance of CirOBl has 
been debat ed in the past with d istance moduli ranging from 
11 .58 mag dMoffat & Voetl 1 19731) to 12.8 mag dVazauez et alJ 
1995I) . 1 3 .00 mag dLortet et al.l 1 1 987l) . and 13.22 mag dLvngal 



1968) 



WR69 is one of the few WC9d stars for which a useful 
IUE spectrum e xists, since the redden ing for this star is rela- 
tively moderate. Willia ms et al.l J2005) measured differences in 
the radial velocities between different observations, possibly in- 
dicating binarity. 

WR 81 is another WC9 star, but one of the few Galactic ex- 
amples without dust emission. Our line fit is remarkably much 
more consistent than for the dusty WC9d stars. 

WR86 is a visual binary classified as WC7 that has a BO 
companion. The star is another example that helps to demon- 
strate which model parameters are obtained when a composite 
spectrum is fitted with a single-star model. We again find that 
the "diluted emission-line" effect places the results from this 
pseudo-fit in the characteristic domain of the log 7Vlog R t -plane 
(Fig.©. 

WR 88 is classified as WC 9 and does not show dust emis- 
sion. Thus it resembles WR 8 1 . Willia ms et al.l ((2005) attribute a 
couple of emission lines to nitrogen, and conclude that WR88 
is either a WC+WN binary or belongs to a previously unob- 
served transitional WN/WC9 subtype. We fitted the spectrum 
with our WC models, and obtained parameters that agree with 
those of other WC9 star s (Fig. [6]). Given the poor quality of the 
iTorres & Massevl ([T987) spectrum we use, we refrain from test 
calculations with enhanced nitrogen abundances. 

W R 90 shows non-thermal radio emission (Chap man et alJ 
119991) . but otherwise no indications of binarity. Our single-star 
fit and the derived parameters are typical of its WC7 type. 

WR 95 has th e spectral subtype WC9d. According to 
Moffat et al] d!977l) . it belongs to the open cluster Trumpler27 
and has a distance of 2 .1 + 0.2 kpc implying a distance modulus 
of P.M . = 11.61 mag. IThe & Stoked d 19701) and iBakker & Thel 
(1983) obtained lower values of 10.17 mag and 11.09 mag, 
which would lead to very low luminosities. 

WR98 was originall y classified as a single star of 
WN8/WC7 transition type. iGamen & Niemelal d2003l) detected 
08-9 type absorption lines superimposed on its spectrum. From 
the SB2 radial velocity curves, they derived the orbital ele- 
ments. We tried to fit the spectrum with our single-star models. 
According to the WN/WC transition type of the spectrum, the 
best-fitting model is basically a hydrogen-free WN model with 
enhanced carbon (5% by mass). However, some spectral features 
cannot be reproduced by any of our single-star models, most 
likely because the contribution of the O-star companion cannot 
be neglected. We therefore consider Fig. B.34 as a "pseudo fit" 
that is in fact unsuitable for this composite spectrum, and omit 
WR 98 from our single-star analysis. 

WR102, also known as Sand 4, is one of two WO stars an- 
alyzed in this work. Classified as W02, it has the same subtype 
as WR 142. To reproduce its spectrum, we needed to use models 
with an enhanced oxygen abundance and very high wind veloc- 



ities. Moreover, the observed emission line profiles of WR 142 
have a round top, which differs from the Gaussian-like shape 
usually encountered in WR spectra. Our models can only repro- 
duce this round shape when we convolve the synthetic spectrum 
for rotational broadening with a vsini of about 1000 km s~'. 
(For wind spectra, flux convolution is only a rough approx- 
imation to account for rotation.) We note that the WO stars 
are very compact, since their radius is smaller than the solar 
radius. Interestingly, the same effect has also been found for 
the most compact W N star in our Galaxy, the WN2 star WR 2 
dHamann et al J 120061) . This rapidly rotating, bare, and compact 
stellar core is certainly a gamma-ray burst ca ndidate. For the dis - 
tance of WR 102, we adopt 3 ± 1 kpc from lDopita et alJ dj990), 
which is based on the nebula G2.4+1.4 that was originally clas- 
sified as a supernova remnant, but is now considered to be stellar 
eje cta from WR 102 d espite this star not being located in its cen- 
ter. iDrew et ail d2004l) calculated a distance of 4.6 kpc based on 
an IR photometry scaling relative to WR 142. We did not use 
this distance as it does not rely on an independent measurement. 
Using a value of 4.6 kpc would infer an extremely high luminos- 
ity of log L/L G = 6.1. 

WR 104 is well known for its pinwheel nebula. It is a bi- 
nary system consisting of a WC9d and a B0.5V star. The latter 
is visually brighter than the WC star dWilliams et al.lll987l). A 
third, fainter component was resolved with HST dWallace et al.l 
2002), hence the WC emission lines are expected to be "diluted". 
Nevertheless, an acceptable fit can still be achieved for WR 104 
with a single-star model that has quite "typical" WC9 parame- 
ters (cf.Fig.[6]l. The distance ( D.M. = 1 1.0 mag) is adopted from 
iLundstrom & Stenholml dl984l) . who assigned WR 104 as a pos- 
sible member to the association SgrOBl. 

WR111 is a prototypical WC5 s tar, which has been fre- 
quently studied. Grafener et al. (2002) basically used the same 
mode ls as the present study, and therefore obtained similar re- 
sults. iGrafener & Hamannl d2005l) constructed a hydrodynami- 
cally consistent model for WR 111, thus showed for the first time 
that WC winds can be explained in terms of radiation-driven 
mass loss. Their model provided a more consistent fit of the line 
spectrum than the semi-empirical models used in the present pa- 
per. The hydrodynamical model of Grafen er & Hamannl (2005) 
has a much higher stellar temperature T* = 140 kK than our 
present study (89 kK). This mainly reflects the different radial 
structures of the two respective models in the deepest zones of 
the wind - we recall that T» is defined as the effective tempera- 
ture corresponding to the radius R*, where the Rosseland optical 
depth reaches 20 ( cf.Eq[Tj). 

The model of IGrafener & Hamannl d2005l) has nearly the 
same luminosity as we obtain from our empirical fit. The mass- 
loss rate of the hydrodynamically consistent model (logM = 
-5.14) is considerably lower than that of our empirical model 
(-4.67). This difference is mainly due to the mu c h high er clump- 
ing contrast assumed by IGrafener & Hamannl (|2005) (D = 50 
in the outer parts instead of 10 in this work). Such strong 
clumping is probably unrealistic, but was needed to compen- 
sate for the incomplete line opacities and achieve sufficient 
radiative driving. A minor part of the mass-loss rate differ- 
ences is due to the terminal wind velocity, which was slightly 
smaller in the hydrodynamic model (2050 km/s) than the value 
adopt ed in the pres e nt pap er (2398 km/s). For the model of 
Graf ener & Hamannl d2005l) . a much higher stellar temperature 
(T* = 140 kK) was assumed than the 89 kK given in the present 
paper. Different stellar radii R* compensate for the effect of lu- 
minosity. The similarity of the emergent spectra demonstrate 
again the parameter degeneracy for very dense winds discussed 
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in Sect. 15. II The star WR 1 1 1 is assumed to be a me mber of the 
SgrOBl association (Lundstrom & Stenholmll 19841) with a dis- 
tance modulus ofD.M. = 11.0 mag. 

WR113 is a WC8d+Q 8-9 binary system 
dCher enashch uk & Karetnikovl l2003h . which has an excess 
in its 2MASS K-band magnitude most probably caused by dust 
emission. The pseudo-fit of its spectrum obtained from our 
single-star models leads to parameters in the binary domain of 
the log 7Vlog/? t -plane (Fig. [6]), obviously due to the dilution of 
the emission lines. 

WR 1 14, classified as WC5, is a me mber of the SerOBl as- 
sociation ( Lundstro m & Stenholm 1984) with a distance modu- 
lus of D.M. = 11.5 mag. The spectrum was found to have diluted 
emission l ines, and the star wa s therefore listed as a binary can- 
didate in lvan der Huchtl(l200ll) . However, our spectral fit and the 
obtained parameters are normal for a single WC5 sta r. 

WR 117 is a "dusty" WC9d star. IWilliams et afl (120051) did 
not find any evidence of an OB companion, and our single- 
star model also fits most features in the observed spectrum with 
typical parameters. The lines are significantly stronger than in 
other WC9 stars. The best-fitting model has a stellar tempera- 
ture of T* = 56 kK, which is relatively high for the WC9 sub- 
class. The position of WR 1 17 in the log 7Vlog 7? t -plane, as well 
as the high terminal wind velocity, are close to or maybe al- 
ready in the WC8 pa rameter region. Some previous papers (e.g. 
IConti & Vaccal [l990l) have indeed classified WR117 as type 
WC8. Unfortunately, our available spectrum does not comprise 
any C n-lines, which would provide the criterion to distinguish 
the subtype WC8 from WC9. 

WR 121 i s another "du sty" WC9d star that is apparently sin- 
gle. IwHliams^eral] (|2005|) could not find any evidence of an OB 
companion, and our single-star model fits most features in the 
observed spectrum with typical parameters. 

WR125 is an SB2 binary system (WC7ed+09III, 
IWilliams et al J [19941) . where the letters "ed" stand for episodic 
dust formation. Our pseudo fit with a single-star model gives 
parameters that are characteristic of composite spectra (cf. 
Fig.©. 

WR 126 shows a unique spectrum that differs from those of 
all other Wolf-Rayet subtypes. Its designation as WC5/WN in- 
dicates that it has predominantly a WC-type spectrum (albeit the 
emission lines are unusually weak), but also relatively strong ni- 
trogen. We tentatively fit the spectrum with a WC-type model 
from our low-carbon (20% mass fraction) grid (Fig. B.46). The 
N iv lines at 7005-703 1 A reveal obviously a significant abun- 
dance of nitrogen, which is not included in our WC models. The 
mass-loss rate is much lower than those typically found for WC5 
stars. A thorough spectral analysis of this transitio n-type star is 
beyon d the scope of the present paper. Following iRadoslavoval 
(1989), the star is probably a member of the Vul OB2 associa- 
tion, which has a distance modulus of D.M. = 13.2 mag. 

WR 132 is classified as WC6 and might be ass ociated with 
an Hi-bubble. Based on that bubble, Arnal| (1 19921) estimated a 
kinematic distance of 4.3 kpc, which we adopt for our luminosity 
scaling. 

WR 135 is one of only four WC8 stars analyzed in this work. 
The star is a possible member of Cyg OB3, lead ing to a distance 
of D.M. = 1 1 .2 mag dGarmanv & Stencel 19921) . The earlier cal- 
culations of ILundstrom & Stenhol m (1984) obtained a slightly 
higher value of 11.6 mag. WR 135 is the only WC8 single star 
with an independent distance estimate. 

WR 137 is a binary system (WC8pd+09) showing periodic 
dust ("pd") formati on. Under the assumpti on that the O star is a 
main-sequence star, IWilliams et al.l d200ll) estimated a distance 



modulus of D.M. = 11.1 mag. The pseudo-fit with a single-star 
model yields parameters that are typical for composite spectra 
(cf. Fig.©. 

WR 142, classified as W02, is the other WO star analyzed 
in this work. It has the same subtype as WR 102. The spectra 
of both stars are very similar, hence so are the results of their 
analyses. We achieved the best fit with a very high terminal 
wind velocity (poo = 5000 km s _1 ) and a very high stellar tem- 
perature of T* = 200 kK. As for WR 102, we had to convolve 
the emergent spectrum with a high rotational broadening veloc- 
ity of 1000 km s' 1 in order to r e produ ce the round shape of the 
emission lines. lOskinova et alJ J2009) detected weak but hard 
X-rays from this object. These X-rays cannot be attributed to 
colliding stellar winds since there is apparently no companion. 
WR 142 is a member of the Berkeley 87 clust er. The distance 
of this cluster has been disputed in the past. iTurner & Forbes! 
(1982) derived a clus ter distance of 946 ±26 pc, which would im- 
ply log L/L Q = 5.5. iMassev et alJ d200ll) and lKnodlseder et all 
(2002) obtained 1.58 kpc and 1.8 kpc, raising the luminosity to 
5.9 or even 6. for the latte r distan ce. In our work, we used the 
latest value of ITurner et alJ d2006l) of a distance of 1230+40pc 
leading to log L/L Q = 5.7. 

WR 143 was reveal ed as a binary (WC4+Be) by 
IVarricatt & Ashokl d2006h . The analysis with our single- 
star models yields parameters that closely fit to the WC 
sequence, possibly indicating that the contribution of the 
Be-type companion to the composite spectrum is relatively 
weak. 

WR144 is one of two Galactic WC4 stars that are not 
suspec ted to be binaries. According to Lundstrom & Stenholm! 
(1984), the star is a possible member of the CygOB2 associa- 
tion with a distance modulus of 11.3 mag, which we adopt for 
our luminosity scaling. This distance implies a luminosity of 
log L/L Q = 5.22. An alternativ e distance modulus of D.M. = 
9.8 mag for Cyg OB2 claimed bv lLinder et al.l d2009l) would lead 
to an implausibly low luminosity of log L/L Q = 4.6. 

WR145 is of the transition type WN7/WC. Its spec- 
trum is nicely fitted by a WN-type model with enhanced 
carbon. The distance modul us of 11.3 mag is taken from 
ILundstrom & Stenholm! d 19841) based on the possible member- 
ship of WR 145 in the Cyg OB 2 association. The star is located 
in a nebula dMiller & Chulll993l) . 

WR146 is a visual bina ry system consis t ing o f a WC5 
star with an 08 companion. Doughe rty et alJ d 19961) also re- 
solved two components with a high-spatial-resolution radio ob- 
servation and classified WR 146 as a colliding wind binary 
due to its non-thermal radio emission. Dougherty et al. (2000) 
performed an in-depth analysis of the different types of radio 
emission from both components and their colliding wind re- 
gion. They suggested that the companion might actually not 
be a single 08 star, but a composite system itself, possibly 
consisting of an 08 and another WC star to explain the high 
mass-loss rate. The poor pseudo fit that can be achieved with 
single-star models yields parameters which are typical of bi- 
naries if the d.e.l. effect is neglected (cf. Fig.©. iDessart et alJ 
(2000) tried to account for the companion's continuum when 
they analyzed this WR star, but still arrived at parameters that 
are atypical of a WC5 star. WR 146 is also listed as a possi- 
ble member of the Cyg OB2 association with D.M. = 11.3 mag 
dLund strom & Stenholm 1 19841). although this membership has 
since been questioned. iDoughertv et all d 19961) obtained D.M. = 
10.4 mag, which would place WR 146 in front of CygOB2 and 
is the value that we used in our pseudo fit. 
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WR154 is a WC6 single star and a possible member of 
the Cep OB 1 association, for which Garma nv & Stencel (1 19921) 
give a distance modulus o f 12.2 mag. Earlier calculations ob- 
tained P.M. = 13.7 mag dSmith et al.l Il990h and 12.53 mag 
dConti& Vaccalll990h . 



